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Abstract 


Saturation  spectroscopy  techniques  were  used  to  demonstrate  the  ability  to  make 
spatially  resolved  temperature  measurements  by  exploiting  the  B3  n(0+u)  -X1  Z(0+g) 
electronic  transition  in  molecular  iodine  (I2).  Doppler-free  saturation  spectroscopy 
measurements  resulted  in  hyperfine  spectral  profiles  for  the  P(70)  17-1  and  P(53)  19-2 
ro-vibrational  transitions  between  the  electronic  levels.  Doppler-limited  saturation 
spectroscopy,  an  adaptation  of  the  Doppler-free  technique  in  which  beams  propagate  in 
the  same  direction  and  cross  within  the  gas  sample,  allowed  for  the  measurement  of 
spatially  resolved  Doppler  profiles  for  the  same  transitions. 

The  profiles  were  measured  at  several  spatial  positions  in  I2  cells  with  measured 
temperature  gradients.  The  technique  did  not  produce  accurate  temperature  readings; 
however,  the  saturation  spectroscopy  variant  did  produce  profiles  that  were  spatially 
resolved,  with  an  interaction  volume  of  approximately  12x1x1  (mm  units). 

This  thesis  was  motivated  by  the  Air  Force's  Airborne  Laser  Laboratory  (ABL), 
which  has  a  need  to  discriminate  a  spatially  distributed  temperature  distribution.  The 
Chemical  Oxygen  Iodine  Laser  (COIL),  which  is  the  high  power  directed  energy  source 
for  ABL,  has  a  low  concentration  of  I2  as  a  reaction  byproduct.  Line-integrated  optical 
temperature  measurements  inside  the  COIL  reaction  chamber  have  been  higher  than 
predicted,  which  can  have  a  dramatic  effect  on  the  threshold  condition  for  the  laser. 
Spatially  resolved  temperature  measurements  would  provide  a  diagnostic  tool  to  identify 
hot  spots  and  heat  gradients  in  the  supersonic  flow  region  of  the  laser. 
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Spatially  Resolved  Temperature  Determination  in  I2  Gas  Using  Doppler- 


Limited  Saturation  Spectroscopy 

I.  Introduction 


A.  Background/History 

Using  absorption  spectroscopy  to  measure  the  temperature  in  a  Doppler-broadened 
gas  sample  is  hardly  a  revolutionary  idea.  In  fact,  the  average  temperature  in  such  a 
medium  can  be  easily  obtained  by  projecting  a  tunable  laser  source  through  the  gas  and 
onto  a  photo-detector.  The  photo-detector  measures  changes  in  absorption  levels  as  the 
source  scans  across  the  entire  base  width  of  a  spectral  line.  The  profile,  assuming  the 
medium  is  primarily  Doppler-broadened,  will  have  a  width  that  is  determined  by  the 
thermal  velocities  of  the  gas  particles  that  interact  with  the  laser  beam.1 

Using  this  method,  the  laser  beam  will  interact  with  all  the  gas  particles  in  the  beam 
path,  and  the  temperature  information  can  be  extracted  from  the  Doppler  width  (Full 
Width  at  Half  Maximum,  FWHM)  using  the  formula:2 


where 


(1) 


Awd  =  Doppler  width  (FWHM) 
c  =  speed  of  light 
T  =  temperature 


coio  =  center  frequency  (rad/s) 
k  =  Boltzmann's  constant 
m  =  molecular  mass 


1 


Equation  (1)  assumes  a  Gaussian  lineshape,  which  is  a  good  assumption  as  long  as 
Doppler  broadening  is  the  primary  broadening  mechanism  in  the  sample.1 

The  major  problem  with  this  approach  for  measuring  the  temperature  is  that  the 
resultant  value  is  the  line-integrated  temperature  over  the  entire  path  of  the  laser  beam 
through  the  sample.  The  measurement  as  such  gives  no  information  about  the  local 
temperature  distribution  in  the  gas  (e.g.  the  measurement  is  not  spatially  resolved),  but 
rather  gives  the  average  temperature  of  the  gas  over  the  laser  beam  path  length. 

In  highly  dynamic  or  turbulent  gases  spatially  resolved  temperature  measurements 
would  provide  useful  information  about  the  media  or  the  flow  dynamics.  Point-to-point 
measurements  of  temperature  in  such  systems  may  differ  significantly  from  the  line- 
integrated  average,  especially  in  situations  where  sharp  temperature  gradients  are 
expected.  Using  a  spatially  resolved  measurement  technique,  a  detailed  mapping  of 
temperature  profiles  in  space  could  provide  diagnosis  of  the  thermal  distribution  of  the 
sample  gas. 

It  is  necessary  that  such  a  temperature  measurement  technique  be  non-intrusive  to  the 
experimental  system.  Of  course,  the  easiest  way  to  achieve  spatially  resolved 

temperature  measurements  is  to  introduce  a  thermocouple  probe  into  the  system.  The 
probe,  if  moved  from  point  to  point,  would  measure  the  resolved  temperatures;  however, 
the  presence  of  the  thermocouple  will  influence  dynamics  and  equilibrium  temperatures 
in  the  gas  species  being  measured.  Optical  measurements  thus  present  an  ideal  means  to 
avoid  significantly  disturbing  the  gas,  while  retaining  the  spatial  resolution. 

The  need  for  the  development  of  a  simple  optical  technique  to  provide  spatially 
resolved  temperature  information  was  inspired  by  the  Airborne  Laser  Laboratory  (ABL) 


2 


technical  development  team,  located  at  Kirtland  Air  Force  Base  in  Albuquerque,  New 
Mexico.  The  ABL  is  using  a  Chemical  Oxygen-Iodine  Laser  (COIL)  for  their  primary 
high  power  laser  source,  and  sees  some  utility  in  being  able  to  measure  a  detailed  spatial 
temperature  map  of  the  gas  inside  the  reaction  chamber  of  the  laser. 

COIL,  being  a  chemical  laser,  relies  on  a  chemical  reaction  to  produce  the  pumping 
energy  required  to  build  a  population  inversion  in  atomic  iodine.  The  primary  reaction 
produces  a  singlet  oxygen  state  (02  a  !A)  with  a  radiative  lifetime  of  about  64  minutes. 
This  state  lifetime  is  not  significantly  diminished  by  collisional  relaxation,  and  therefore 
acts  as  a  metastable  energy  reservoir.  Atomic  iodine,  which  dissociates  from  I2 
downstream  from  this  reaction,  has  an  I  2Pi/2  energy  state  that  is  close  in  energy  to  the 
singlet,  excited  oxygen  state.  The  energy  from  the  metastable  reservoir  state  of  the 
oxygen  can  couple  to  excite  the  two-level  inversion  in  the  atomic  iodine.3  Figure  1 
shows  a  conceptual  comparison  of  energy  levels  for  the  primary  components  in  the  COIL 
reaction  chamber. 

The  resonant  energy  transfer  reaction  is  given  by: 4 


o2('a)  +  I(2P3/2)  ®  02(3Z)  +  I'(2Pi/2). 


(2) 


Figure  1.  Energy  level  comparison  for  I,  02,  and  I2.  The  curved  arrow  shows  the 
coupling  between  the  excited  oxygen  state  and  the  atomic  iodine  upper  lasing  state. 
The  primary  lasing  transition  in  I  is  also  designated,  with  a  wavelength  of  1.315  pm. 
Energy  levels  of  I2  are  included  to  show  the  B-X  transitions  investigated  in  this 
thesis.3 

The  reaction  cell  for  the  COIL  employs  a  vacuum  to  force  the  oxygen  through 
supersonic  nozzles.  Molecular  iodine  (I2)  is  injected  at  the  entrance  plane  of  the 
supersonic  nozzles,  as  shown  in  Figure  2.  The  nozzles  are  present  to  cool  the  gases  as 
they  travel  through  the  laser  cavity,  in  order  to  control  temperature-dependent  losses  and 
to  increase  the  power  flux.  The  I2  dissociates  in  a  complex  series  of  kinetic  interactions 
with  the  singlet  oxygen  state;  however,  trace  concentrations  of  I2  still  remain  past  the  exit 
plane  of  the  nozzles.3  The  I2  will  likely  have  the  same  temperature  as  the  surrounding 
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gas  components  as  the  flow  reaches  the  resonant  cavity;  it  is  this  downstream  temperature 
which  is  of  interest. 


Laser  Output 


Figure  2.  COIL  nozzles  and  optical  cavity.  The  L  is  injected  through  the 
supersonic  nozzles  and  dissociates,  allowing  for  the  resonant  energy  transfer  from 
the  excited  oxygen  to  the  atomic  iodine.  The  lasing  transition  occurs  aft  of  the 
nozzles  in  the  resonant  laser  cavity.3 

The  temperature  of  the  gas  at  the  exit  plane  of  the  nozzles  determines  the  threshold 
condition  for  lasing;  a  temperature  change  from  room  temperature  to  160  K  reduces  the 
ratio  of  singlet  oxygen  to  total  oxygen  concentration  required  for  lasing  from  15%  to 
5%.3  The  efficiency  is  therefore  greatly  reduced  by  higher  temperatures;  since  the  COIL 
is  to  be  the  primary  laser  for  development  in  ABL,  efficiency  is  paramount  in  system 
design  priority  for  the  development  of  a  practical  platform.  Weight  of  fuel  gases,  heat 
production,  and  laser  output  power  are  all  contingent  upon  the  efficiency  of  the  system 

Informal  communications  between  the  Air  Force  Institute  of  Technology  (AFIT)  and 
the  Air  Force  Research  Laboratory’s  ABL  branch  (AFRL/DEL)  have  indicated  that  line- 
integrated  temperature  measurements  made  of  the  I2  trace  gas  at  the  nozzle  exit  point 
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have  yielded  unexpected  results.  The  measured  temperature  within  the  reaction  chamber 
of  the  laser  was  measured  50  degrees  (Kelvin)  higher  than  predicted. 

The  ability  to  resolve  and  map  the  spatial  temperature  variation  in  the  COIL  will 
assist  the  ABL  in  determining  the  source  of  the  high  line-integrated  temperature 
measurement.  It  is  even  possible  that  there  is  not  a  problem  with  the  temperature 
distribution  at  the  nozzle  exit  point,  but  that  the  higher  measurement  is  a  result  of  the  line 
integration  in  the  interaction  volume.  In  any  case,  it  is  necessary  to  develop  a  technique 
that  will  assist  in  the  diagnosis  of  the  problem,  or  to  determine  if  a  problem  even  exists. 

After  several  discussions  with  AFRL/DEL,  AFIT  personnel  suggested  that  a  variant 
of  a  well-known  nonlinear  laser  spectroscopy  technique  could  provide  spatially  resolved 
temperature  information.  This  technique,  known  as  saturation  spectroscopy,  is  generally 
applied  when  one  wishes  to  measure  the  hyperfine  spectrum  that  comprises  a  Doppler- 
broadened  transition.  Counterpropagation  of  two  split  beams  from  a  tunable  source 
within  the  gas  sample  allows  the  elimination  of  the  Doppler  broadening  altogether  via 
Doppler-free  saturation  spectroscopy. 

The  elimination  of  the  primary  broadening  mechanism  allows  the  component 
hyperfine  line  positions  and  widths  to  be  measured.  Figure  3  shows  the  comparison  of 
the  absorption  line  to  its  component  hyperfine  line  structure  for  the  P(70)  17-1  transition 
in  h.  The  broadened  profile  is  in  actuality  several  orders  of  magnitude  greater  in 
amplitude  than  the  hyperfine  lines.  The  hyperfine  line  intensities  have  been  exaggerated 
for  the  sake  of  comparison.  The  additions  of  these  15  hyperfine  lines,  when  Doppler 
broadened,  make  the  smooth  profile. 
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Wave  Number  (cm-1) 


Figure  3.  P(70)  17-1  Absorption  line  with  component  hyperfine  structure.  There 
are  15  hyperfine  lines  for  this  transition  that,  when  Doppler  broadened,  sum 
together  to  form  the  broadened  transition.  This  diagram  is  a  relative  illustration, 
since  the  broadened  transition  is  in  actuality  several  orders  of  magnitude  greater  in 
magnitude  than  the  hyperfine  lines. 


In  the  Doppler-limited  modification  to  the  Doppler-free  technique,  the  two  laser 
beams  instead  copropagate  and  cross  within  the  sample  to  yield  the  Doppler  broadened 
profile  only  at  the  point  where  the  two  beams  cross.  If  the  beams  cross  at  a  small  angle, 
there  is  again  interaction  with  the  same  velocity  group  and  the  temperature  information 
that  is  extracted  from  the  Doppler  width  of  the  resulting  line  is  spatially  resolved.  The 
position  at  which  the  beams  cross  can  also  be  altered  to  obtain  the  temperature 
dependent,  spatially  resolved  profiles  at  different  spatial  locations  in  the  sample. 
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B.  Objective 


The  primary  objective  of  this  work  was  to  demonstrate  the  capability  of  Doppler- 
limited  saturation  spectroscopy  for  making  spatially  resolved  temperature  measurements 
in  a  molecular  iodine  gas.  This  included  obtaining  an  understanding  of  the  hyperfine 
interactions  that  compose  I2  absorption  lines,  and  applying  principles  of  Doppler 
broadening  to  retrieve  temperature  information  from  the  profiles  of  the  measured  lines. 
The  Doppler-free  saturation  spectroscopy  experiment  was  slightly  modified  to  not 
exclude,  but  to  exploit  the  spatially  resolved  Doppler  information  to  extract  temperatures. 

The  laser  source  for  the  experiment  was  the  Coherent  899  Ring  Dye  Laser,  with  a 
Rhodamine  590  dye  to  allow  for  tuneability  from  the  570  to  6 1 0  nm  range.  The  overall 
approach  was  to  measure  the  hyperfine  structure  for  several  broadened  transitions,  which 
were  chosen  because  of  their  isolation  in  frequency  space  from  other  absorption  lines, 
and  for  characteristics  that  may  potentially  increase  sensitivity  to  temperature  variations. 
Hyperfine  measurements  were  made  using  standard  Doppler-free  saturation  spectroscopy 
techniques.  The  experimental  setup  was  then  modified  for  Doppler-limited  spatially 
resolved  measurements  by  changing  the  direction  of  propagation  for  the  stronger  of  two 
laser  beams  (the  pump  beam). 

The  resulting  Doppler-limited  profiles  were  then  fit  starting  with  the  component 
hyperfine  positions  and  natural  linewidths  as  measured.  A  Doppler  width  was  then 
assigned  to  these  lines  and  was  allowed  to  change  as  one  of  the  fit  parameters;  therefore 
the  broadened  profile  could  be  recreated  by  the  addition  of  the  hyperfine  lines  with  a 
shared  Doppler  width.  The  profiles  were  measured  at  different  points  along  the  length  of 


8 


a  cell,  heated  at  one  end  in  order  to  create  a  measurable  temperature  gradient  in  the  gas 
within  the  cell.  These  measurements  were  compared  to  those  made  by  introducing  a 
thermocouple  probe  into  the  cell  under  identical  pressure  conditions. 

C.  Presentation 

The  relevant  theoretical  background  information  is  presented  in  Chapter  II.  The 
development  covers  the  electronic  transitions  in  molecular  iodine,  hyperfine  splitting, 
Doppler-broadened  transitions,  Doppler-free,  and  Doppler-limited  saturation 
spectroscopy.  Chapter  III  describes  experimental  assemblies  and  includes  the  procedures 
for  hyperfine  measurements  and  for  Doppler-limited  measurements  using  two  Iodine  cell 
designs.  Chapter  IV  presents  the  experimental  results,  including  hyperfine  structure 
profiles  and  corresponding  Doppler  broadened  profiles,  which  are  presented  with  detailed 
fit  parameters  and  fit  residual  analysis.  Chapter  V  gives  the  discussion  of  the  results  and 
recommendations,  which  include  an  outline  of  possible  future  work  in  this  area  of 
interest. 
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II.  Theory 


A.  Electronic  Transitions  in  I2 

The  absorption  spectrum  of  molecular  iodine  has  been  well  investigated  and 
catalogued  since  the  birth  of  molecular  spectroscopy.  Specifically,  the  ro-vibronic 
transitions  between  the  X1I+g  B3n+U  electronic  states  have  been  extensively 
catalogued  by  using  high  resolution  Fourier  transform  spectroscopy,  and  are  presented  in 
an  atlas  by  Gerstenkom  and  Luc  to  an  accuracy  of -0.002  cm'1.4  A  portion  of  the  I2 
spectrum  is  shown  in  Figure  4;  the  primary  lines  chosen  for  the  experiment  lie  within  this 
range.  The  wealth  of  empirical  information  concerning  absorption  transitions  in  I2  have 
led  to  a  fair  understanding  of  the  physical  stmcture  of  the  molecule,  and  of  its  energy 
level  structure. 


Figure  4.  Iodine  absorption  spectrum.  This  spectrum  shows  the  relative 
transmission  through  12  as  a  monochromatic  laser  source  is  scanned  across 
frequency  space.5 
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As  with  any  diatomic  molecule,  molecular  iodine  has  rotational  and  vibrational 
degrees  of  freedom.  Rotational  and  vibrational  states  are  created  by  the  quantized  nature 
of  these  motions,  essentially  splitting  the  electronic  state  as  the  physical  configuration  of 
the  molecule  changes.8  An  illustrative  potential  diagram,  Figure  5,  shows  the  typical 
energy  state  hierarchy  for  a  diatomic  molecule.  Within  each  electronic  state  there  are 
vibrational  states,  in  which  there  are  rotational  states.  Simple  transitions  from  one 
rotational  state  to  another  within  the  same  vibrational  level  are  possible;  furthermore, 
vibrational  transitions  are  possible  within  the  same  electronic  state.  Electronic 
transitions,  however,  are  generally  accompanied  by  ro-vibrational  transitions.7 

Vibration  in  the  molecule  causes  the  internuclear  separation  to  change 
anharmonically.  The  motion  is  not  perfectly  harmonic  because  the  Coulomb  interaction 
between  nuclei  is  stronger  as  the  internuclear  separation  decreases,  and  is  weaker  as  the 
separation  increases.  This  causes  the  electronic  potential  curves  to  be  asymmetric  when 
compared  to  the  theoretical  parabolic  potential  of  the  harmonic  oscillator.11  Vibrational 
states  are  determined  by  the  shape  of  the  anharmonic  potential,  commonly  referred  to  as 
the  Morse  potential.12 

The  term  value  for  the  vibrational  energy  levels  is  given  by  the  Dunham  Expansion, 
which  is:11 

Jmax 

Gv=  Z  Yo/ 

j  =  0 

where  v  is  the  vibrational  quantum  number  and  the  expansion  coefficients  (Y0i,Y02,etc.) 
are  empirically  determined  by  robust  fitting  routines  that  examine  numerous  vibrational 
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transitions  in  the  molecule.  Yoo  is  the  electronic  energy  level,  determined  at  the 
equilibrium  point  on  the  Morse  potential  curve;  this  expansion  coefficient  is  often 
denoted  as  Te.  In  iodine,  these  coefficients  are  very  well  determined  and  are  used  to 
predict  and  classify  transitions.  Table  la  shows  the  expansion  coefficients  for  the  ground 
state  (X1  l(0+g)),  while  those  for  the  B3  flf0+u)  state  are  shown  in  Table  lb. 


Figure  5.  Typical  diatomic  energy  level  structure.  Each  electronic  level  has 
vibrational  levels,  and  each  vibrational  level  has  rotational  levels.9 


12 


For  simple  rotation,  considering  the  molecule  to  be  a  rigid  rotator,  one  finds  that  the 
resulting  energy  levels  are  evenly  spaced  as  defined  by  the  term  symbol: 

F(J)=B  J(J+1)  (4) 

where  B  is  the  rotational  constant  and  J  is  the  rotational  quantum  number.  In  actuality, 
since  the  molecule  is  not  completely  rigid,  there  is  a  centrifugal  distortion  that  must  be 
taken  into  account.  The  term  value  is  then: 

F(J)=B  J(J+1)+D  J2(J+1)2+H  J3(J+1)3+...  (5). 

Rotational  branch  structure  is  formed  as  a  consequence  of  the  selection  rules  for 
transitions  in  diatomic  molecules.  Transitions  from  state  to  state  can  occur  in  a  diatomic 
molecule  if  AJ  =  -1,  0,  or  +1;  the  branches  are  defined  as  follows: 

J'-J"  =  OQ-Branch 
T  -  J"  =  +1  R-Branch 
J  -  J"  =  -1  P-Branch 

Figure  6  illustrates  branch  structure  in  typical  diatomics  by  identifying  the  specific 
transitions  and  frequency  distribution  for  a  simple  case. 

The  identification  of  the  rotational  transitions  is  very  important  for  this  work,  since 
the  number  of  hyperfine  lines  in  a  transition  are  determined  by  the  rotational  quantum 
number  in  the  ground  state,  J".  In  I2  there  is  no  Q  branch  formation  because  the  resultant 
angular  momentum  (Q)  for  each  state  is  zero.  This  is  according  to  the  selection  rule  13 


AJ  =/  0  when  O'  =  O"  =0 


(6) 
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As  one  can  expect,  vibration  and  rotation  occur  simultaneously  in  a  molecule.  If  the 
molecule  vibrates,  the  intemuclear  distance  is  changing;  therefore,  the  angular 
momentum  of  the  molecule  is  changing.  The  reverse  is  true  as  well.7  In  order  to 
compensate  for  this  mutual  interaction,  the  rotational  constants  are  themselves 
expansions  that  depend  on  powers  of  (v+1/2).14 


P  branch'— — — |  j— —  R  branch 
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Figure  6.  Branch  structure  in  electronic  transitions.  The  P-and  R-branch 
transitions  are  conceptually  illustrated  for  rotational  transitions  from  the  excited 
(Single  -  primed)  electronic  state  to  the  lower  (Double  -  primed)  state.11 
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The  possible  energy  levels  (in  cm'1)  can  then  be  approximated  by  the  set  of  equations 
labeled  (7)  below:15 


A  sufficient  set  of  expansion  coefficients  have  been  determined  for  the  B-X 


transition  in  molecular  iodine  and  are  presented  in  Table  1,  and  are  identified  by  the  i  and 
j  subscripts  corresponding  to  Equation  (7).  Since  the  lower  vibrational  levels  of  the 
ground  state  are  more  densely  populated  states  at  room  temperature,  the  ground  state 
parameters  are  only  accurate  (within  0.01  cm'1)  for  the  first  10  vibrational  levels.14 
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Table  la.  Expansion  parameters  for  the  X1  S(0+g)  state  of  I2 .  Valid  for  v"=0-914 


— 

1— 

iH 

ijlli 

6 

3.4368670E-02 

4.56812301-09 

2.1451860E+02 

-1.1398970E-04 

2.1204390E-11 

-6.0722840E-01 

-2.72029 10E-07 

6.7374060E-13 

-1.3861 100E+01 

-5.4477080E-09 

Table  lb.  Expansion  parameters  for  the  B3  n(0+u)  state  of  I214 


— 

— 

QB|B| 
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iBii 

15769.0485  (Te) 

2.8993694599E-02 

6.125767E-09 

-2.150047340E-15 

■■I# 

125.6724 

-1.406799398E-04 

1.418420E-10 

7.915679522E-15 

.  .m . 

-7.5267700E-01 

-5.088972976E-06 

-2.825041E-12 

-1.023966180E-14 

■. . ML . 

-3.2462820E-03 

8.751 145403E-07 

5.544495E-13 

4.852619270E-15 

:  (Mi  - 

1.8757360E-05 

-1.171736403E-07 

-2.644252E-14 

1.19975 1560E- 15 

-3.4141240E-06 

9.702596783E-09 

9.539521E-16 

1.792143 124E-16 

■  « 

2.0049980E-07 

-5.304859892E-10 

-1.864100E-17 

-1.755753490E-17 

— 

-6.9504140E-09 

1.9571 12328E-11 

1.813907E-19 

1.185081818E-18 

1.5178990E-10 

-4. 89486405 0E- 13 

-6.053203E-22 

-5.675954430E-20 

lltlll  {$*$} 

-1.8997900E-12 

8.1736417235E-15 

1.958960707E-21 

WmmmMjmiM 

1. 22694 10E- 14 

-8.715621632E-17 

-4.888492740E-23 

-3. 1373 120E-17 

5.360997383E-19 

8.739722535E-25 

-1.445945450E-21 

- 1.091 263 060E-26 

9.03491 1061E-29 

-4.455602540E-31 

9.9040298 15E-34 

Once  the  parameters  are  used  in  their  corresponding  expansion  terms,  the  energy 


level  structure  can  be  approximated.  The  first  few  energy  levels  in  I2  are  shown  in  Figure 
7.  There  is  a  dissociative  (repulsive)  energy  level  which  intersects  the  B3  IT(0+u)  state 


which  causes  significant  predissociation.  There  are  many  strong  transitions  in  molecular 
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iodine  in  the  visible  spectrum,  caused  by  the  strong  absorption  character  created  by  the 
difference  of  the  equilibrium  separation  in  the  X  and  B  bound  states.16 


the  B  state  has  significant  predissociation  caused  by  the  crossing  of  a  repulsive 
potential  level. 16 

The  importance  of  identifying  specific  transitions  is  made  clear  in  the  next  section. 


The  transitions  for  this  experiment  were  largely  chosen  because  the  absorption  lines  were 
well  isolated  in  frequency  space.  Once  chosen,  line  assignments  had  to  be  made  using 
the  expansion  parameters  and  the  preceding  development  in  order  to  determine  the 
primary  characteristics  of  the  hyperfine  structure  for  each  absorption  line. 
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B.  Hyperfine  Structure  in  I2 


Prominent  absorption  lines  in  a  low  pressure  molecular  iodine  gas  sometimes  have 
FWHM's  that  are  as  much  as  twice  as  large  as  that  calculated  for  a  simple  Gaussian  line 
shape  determined  by  the  Doppler  effect.17  The  discrepancy  had,  prior  to  1970,  been 
attributed  to  unresolved  hyperfine  splittings.  The  advent  of  Doppler-free  saturation 
spectroscopy  allowed  for  the  resolution  of  hyperfine  interactions,  which  previously  had 
been  measured  and  characterized  only  for  transitions  in  the  microwave  regions  of  the 

17  18 

spectrum.  ’ 

Hyperfine  splittings  are  attributed  to  the  interaction  between  the  electronic  and 
nuclear  spin  angular  momentum,  which  splits  the  fine  spectrum  created  by  spin-orbit 
coupling.  The  nucleus  of  I127  has  a  large  magnetic  dipole  moment,  caused  by  the  spin  of 
the  nuclear  particles;  in  addition,  there  is  a  large  electric  quadrupole  moment  created  by 
the  nonuniform  charge  distribution  in  the  nucleus.  Like  all  spectroscopic  transitions, 
hyperfine  splittings  are  caused  by  electric  and  magnetic  fields  that  interact  with  these 
moments. 10 

The  quantum  number  J  is  typically  defined  by  the  addition  of  the  spin  and  orbital 
angular  momentum.  If  we  denote  the  resultant  spin  due  to  the  charged  particles  within 
the  nuclei  in  a  diatomic  particle  as  T,  then  we  can  define  a  new  quantum  number  for  the 
total  angular  momentum,  inclusive  of  nuclear  spin:13 


F=J  +  T,J  +  T-1,J  +  T-2,...|J-T|  (9) 

T=2I,  where  I  is  the  spin  vector  of  each  nucleus 
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I127  has  a  nuclear  spin  of  5/2,  such  that  the  possible  values  for  T  are  5,  4,  3,  2,  1,  0.17 
The  number  of  hyperfine  states  corresponding  to  each  rotational  state  is  determined  by 


the  sum  of  possible  T  degeneracies,  and  is  mathematically  given  by: 

NsV  2-Tfl 

s  La  n 

n 

(10) 

There  is  a  restriction,  however,  on  the  possibilities  for  T-values,  depending  on 
whether  the  states  have  an  even  or  an  odd  J.19  The  total  wavefunction  of  a  homonuclear 
diatomic  molecule  with  any  half-integer  nuclear  spin  must  be  antisymmetric  under 
exchange  of  the  nuclei,  in  agreement  with  the  Pauli  exclusion  principle.  The  B  state  in 
molecular  iodine  is  ungerade,  meaning  that  it  has  odd  electronic  symmetry,  while  the  X 
state  is  gerade.20 

In  an  ungerade  state,  such  as  the  B  state  in  h,  only  states  of  odd  T  values  can  combine 
with  odd  J  states;  therefore,  odd  J-states  have  Ns=21  by 

[2  (5)  +  1]  +  [2  (3)  +1]  +  [2  (1)  +1]=1 1+7  +  3=21 

Likewise,  in  an  ungerade  state,  the  even  levels  can  only  combine  with  even  T  values, 
such  that  NS=T5.  Because  of  the  change  in  electronic  symmetry  for  gerade  states,  such  as 
the  X  state,  the  rules  are  reversed.  For  odd  J  levels  in  gerade  states  there  are  21  substates 
and  for  even  levels  there  are  15. 10 

Therefore,  for  P-  or  R-branch  transitions,  there  are  the  same  numbers  of  splitting 
states  in  both  of  the  energy  levels,  such  that  there  are  either  21  or  15  hyperfine  lines  in 
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one  broadened  absorption  line.  For  even  J",  there  are  15  hyperfine  lines,  while  for  odd 
J",  there  are  21  hyperfine  lines  for  I2 .  The  selection  rule  for  hyperfine  transitions  is  that 
AF=0,  meaning  that  there  will  be  as  many  transitions  as  there  are  splittings  for  a  given  J'- 
J"  transition  for  I2.11  A  conceptual  level  diagram  is  shown  in  Figure  8;  the  dotted 
transitions  are  those  excluded  when  there  are  only  15  hyperfine  lines  splittings  in  the 
states. 

There  is  no  electronic  magnetic  field  along  the  J  direction  in  either  of  the  two  states 
being  investigated.  The  B  state  may  have  a  magnetic  moment  which  is  along  the 
internuclear  axis;  however,  since  Q  =  0,  the  direction  of  this  field  must  be  orthogonal  to 
J.18  Therefore  the  hyperfine  splitting  caused  by  the  magnetic  dipole  moment  interaction 
is  a  minor  effect,  caused  by  the  interaction  of  the  nuclear  spin  with  a  component  of 
angular  momentum  induced  by  the  end  over  end  rotation  of  the  molecule  as  a  whole.21 
The  electric  quadrupole  term  is  the  dominant  mechanism  behind  the  hyperfine  splitting  in 
I2 .  The  total  Hamiltonian  for  the  splitting  terms  is  Hneq  +  Hm,  where  the  terms  are  of 
the  form10 


HNEO=CQqj' 


3<vjM(vj)- 


2T^V  I) 


(11) 


where 

Q=electric  quadrupole  moment  qj=average  E-field  gradient  along  J  direction 
e=charge  of  an  electron 


20 


where  |iG/Ii=Spin-rotation  coupling  constant 


Figure  8.  Hyperfine  transitions  in  I2.  In  each  electronic  state  the  individual 
splittings  are  due  to  the  possible  quantized  F-values.  The  dotted  lines  designate 
transitions  that  do  not  occur  for  even  J"  interactions,  such  that  there  are  21  odd  J 
hyperfine  transitions  and  15  even  J"  hyperfine  transitions.19 
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In  actuality,  finding  the  energy  values  using  the  Hamiltonian  equations  (11)  and  (12) 
requires  diagonalization  of  a  nontrivial  matrix.  However,  an  approximate  solution  can  be 
formulated,  with  the  assumption  that  J  »  I;  this  condition  will  be  met  for  the  transitions 
in  this  experiment.  Recognizing  that  M  is  the  projection  of  I  along  the  axis  parallel  to  F, 
the  total  angular  momentum  vector,  the  approximate  solution  of  the  equations  takes  the 
form:10 


eQdj 

EnEQ”  8-Ij  ^-l!  -  l) 


3-VM!2  +  M22J-hj 


M! 


Mr(M,  +  l)- Ij-^  +  i 


■h 


Mr 


M2-(M2  +  l)-I1-(l1+l)-hi 


(13) 


E  = 


m 


(Mi  +M2)-J 


(14) 


These  equations  provide  approximate  expressions  for  relative  position  of  the 
hypcrfine  lines.  The  theoretical  intensities  of  the  lines,  however,  are  proportional  to 
2F+1.18  It  is  important  to  note  that  the  constant  terms  at  the  beginning  of  each  equation 
are  transition  dependent.  For  this  work  the  major  requirement  is  that  the  hyperfine 
intensities  are  matched  with  the  correct  relative  line  positions;  absolute  line  positions  are 
not  needed  for  the  analysis.  Therefore,  in  this  thesis,  I  will  choose  the  representative 
values  eQqo=-2900  MHz,  pG/Ii=49.0  kHz,  AeQq=1936  MHz.  These  are  the  values  used 
by  Levenson  for  the  P(  117)  21-1  transition;  the  use  of  these  values  will  give  the  relative 
positions  of  the  hyperfine  lines.10 

Figure  9  shows  the  typical  electric  quadrupole  hyperfine  splittings  in  I2  versus 
increasing  J.  Again,  the  dotted  lines  show  the  transitions  which  do  not  occur  for  even  J" 
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transitions,  in  which  there  are  only  15  hyperfine  lines.  The  magnetic  dipole  splitting 
tends  to  shift  these  lines,  and  can  cause  certain  components  to  cross  over  one  another. 
For  the  purpose  of  this  thesis,  both  coupling  terms  will  be  included  in  the  approximate 
solution  to  the  Hamiltonian. 


Figure  9.  Hyperfine  splittings  from  the  electric  quadrupole  interaction  in  I2.  As  J 
increases,  the  splittings  are  closer  to  one  another  in  frequency  space.  Av  is  the 
relative  frequency,  when  an  arbitrary  point  in  the  hyperfine  spectrum  is  chosen  as  a 
zero,  or  center  point.10 


The  addition  of  the  hyperfine  components  forms  a  single  broadened  absorption  line 
and  account  for  the  discrepancy  in  the  FWHM  that  was  mentioned  at  the  beginning  of 
this  section.  The  Doppler  broadening  of  the  hyperfine  lines  for  different  temperatures 
make  profiles  that  are  unique  in  the  sense  that  each  profile  belongs  to  a  separate 
temperature.  Figure  10  illustrates  the  theoretical  superposition  of  the  hyperfine  lines  at 
increasing  temperatures.  The  plot  starts  with  veiy  sharp  hyperfine  lines  at  the  bottom  of 
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the  plot.  Each  subsequent  curve  in  the  vertical  direction  constitutes  a  50-degree  (Kelvin) 
increase  in  temperature;  all  relative  hyperfine  intensities  are  theoretically  justified  by  a 
calculation  of  the  factor  2F+1  for  each  line. 


Figure  10.  Single  transition  profiles  vs.  temperature.  The  hyperfine  intensities  are 
set  to  be  proportional  to  2F+1.  At  0  degrees  Kelvin,  the  hyperfine  lines  form  a 
discrete  spectrum,  as  in  the  bottom  of  this  plot.  As  temperature  increases,  the 
profile  takes  a  less-defined  shape  as  the  broadened  hyperfine  lines  begin  to  overlap. 
Each  subsequent  curve  represents  a  50  Kelvin  increase  in  temperature. 

C.  Doppler  Broadened  Transitions 

The  frequency  ®o  which  interacts  with  a  particle  moving  at  velocity  v  is  Doppler 
shifted  in  the  particle's  frame  of  reference  to  a  new  frequency  cot,  +  k  •  v,  where  k  is  the 
propagation  vector  of  the  light.  Generally,  when  dealing  with  light  interacting  with  a  gas 
of  atoms  or  molecules,  any  volume  will  have  a  Maxwellian  distribution  of  velocities  in 
all  directions. 
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The  absorption  coefficient  for  the  gas  being  illuminated  is  typically  of  the  Voigt 
integral  form:22 


so-h-(2-K-T) 


(15) 


y=transverse  relaxation  rate  y  =  (yl  +  y2)/2 
x=longitudinal  relaxation  time  x  =  (yl'1  +  y2'1)/2 
(i=electric  dipole  operator  s0=permittivity  of  free  space 

K=Boltzmann's  constant  y=the  characteristic  Rabi  frequency 

k=propagation  constant  T=temperature 

The  Voigt  integral  is  essentially  the  convolution  of  a  Gaussian  lineshape  and  a  Lorentzian 
lineshape,  in  this  case  the  Doppler  broadened  profile  and  the  homogeneously  broadened 
profile.  In  the  limiting  case,  when  the  Doppler  broadening  is  much  more  prominent  than 
the  homogeneous  broadening,  the  absorption  coefficient  reduces  to:22 
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and  the  Doppler  width  (FWHM)  is  given  as:23 
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From  a  measurement  of  the  FWHM  of  a  single  broadened  line,  one  can  extract  the 
temperature  information,  simply  from  the  shape  of  the  distribution. 

In  the  case  where  the  absorption  line  is  composed  of  hyperfine  lines,  as  in  this 
experiment,  the  broadened  absorption  profile  has  a  FWHM  that  is  much  larger  than  what 
is  predicted  by  Equation  (1).  This  is  because  the  profile  is  actually  the  sum  of  the 
hyperfine  fines,  each  of  which  is  broadened  according  to  the  temperature  of  the  I2  gas. 
When  trying  to  extract  temperature  information,  it  is  necessary  to  incorporate  the 
positions  of  the  hyperfine  components.  Their  Doppler  widths,  barring  any  nonlinear 
phenomena  in  the  medium,  are  shared  and  reflect  the  true  temperature  of  the  gas. 

If  one  wishes  to  illuminate  the  medium  with  a  monochromatic,  coherent  fight  source: 

E  (r  ,  t)  =  Eo  cos  (co  t  +  n  -  k  •  r)  (17) 

k=propagation  vector  r=3  -D  position  vector 

co=frequency  of  oscillation  t=time 


only  those  molecules,  which  are  Doppler,  shifted  into  resonance  with  the  incident  field 
will  interact.23  The  incident  radiation  can  therefore  only  be  absorbed  by  molecules 
within  a  narrow  region  of  velocity  space.  Molecules  that  fie  outside  this  region  see  a 
frequency  shifted  so  far  from  resonance  that  absorption  is  impossible.10  The  population 
of  the  upper  state  (N2)  will  be  increased  at  the  resonant  velocity  component  vres:22 
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Likewise,  the  lower  state  (Nl)  population  decreases  at  the  same  velocity  component.  The 
velocity  space  profiles  for  the  upper  and  lower  state,  then,  have  distributions  that  appear 
similar  to  Figure  1 1 . 
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Figure  11.  Velocity  distributions  for  Nl(v),  N2(v).  The  incident  radiation  interacts 
with  one  velocity  group,  causing  absorption  to  occur  only  at  vres.  This  causes  Nl  to 
be  depleted  at  this  velocity,  and  N2  to  increase.  22 

This  phenomenon,  known  as  hole  burning,  occurs  because  a  Doppler-broadened 
profile  can  be  described  by  an  addition  of  many  naturally  broadened  lines,  each  with  a 
Doppler-shifted  center  frequency,  as  in  Figure  12. 


Figure  12.  Doppler  profile  as  a  superposition.  The  broadened  line  is  composed  of 
lines  that  are  naturally  broadened,  each  with  a  center  frequency  that  is  Doppler 
shifted  with  respect  to  the  incident  radiation.24 
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The  change  in  the  velocity  distributions  for  N1  and  N2,  illustrated  in  Figure  1 1, 
cause  a  hole  to  form  in  the  absorption  profile  at  the  source's  frequency.  At  this  point, 
called  the  Bennett  Hole,  the  absorption  in  the  medium  is  lower  because  the  upper  states 
become  more  populated.  In  fact,  if  the  source  power  is  pushed  past  the  saturation 
intensity  for  the  medium,  the  medium  becomes  transparent  to  the  source  frequency.  If 
the  source  is  scanned  over  the  frequency  range  of  the  entire  broadened  transition,  the 
Bennett  hole  will  move  through  the  entire  profile  in  frequency  space. 

The  frequency  of  the  monochromatic  source  therefore  defines  what  velocity 
component  in  the  medium  will  interact  with  the  beam  Figure  13  shows  three 
possibilities  for  frequency  interaction  in  a  Doppler  broadened  medium.  For  case  a,  the 
light  can  only  interact  with  the  molecules  that  have  no  projection  of  their  physical 
velocities  onto  the  propagation  vector;  the  particles  that  interact  with  this  incident 
frequency  have  velocities  which  are  perpendicular  to  the  incident  beam  direction.  When 
v  >  v0  (case  b),  the  wave  interacts  with  molecules  that  have  positive  projections  of  then- 
velocities  onto  the  wave  propagation  vector,  and  a  hole  is  burnt  into  one  side  of  the 
absorption  curve.  Likewise,  case  c  involves  a  negative  projection  of  particle  velocity 
onto  the  propagation  vector,  and  the  hole  appears  at  the  opposite  side  of  the  profile  from 
case  b. 
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Figure  13.  Interactions  for  three  incident  frequency  cases.  The  three  distinct 
frequencies  interact  with  molecules  that  have  different  axial  velocities,  causing  holes 
to  be  burnt  at  different  locations  in  the  Doppler-broadened  profile.24 


Now  it  is  useful  for  the  purpose  of  this  thesis  to  consider  the  case  where  two 
monochromatic  waves  of  the  same  frequency  are  allowed  to  counterpropagate  through 
the  medium.  The  beam  in  one  direction  bums  a  hole  on  one  side  of  the  absorption 
profile,  while  the  beam  in  the  other  direction  burns  a  hole  on  the  other  side;  there  will  be 
two  holes  that  are  symmetric  about  the  center  frequency  of  the  absorption  curve.  The 
beams  interact  with  separate  (in  fact,  opposite)  axial  velocity  groups,  and  so  there  are  two 
distinct  depressions  in  the  absorption  profile.23 

If,  however,  the  source  is  scanned  across  the  profile,  the  beams  will  reach  a  point 
where  they  both  cross  the  center  frequency.  At  this  point  the  beams  are  interacting  with 
the  same  velocity  subgroup,  namely  vz=0,  as  Figure  14  a  illustrates. 


Figure  14  a.  The  Bennett  hole  caused  by  counterpropagation.  Two  beams  of  the 
same  frequency  are  allowed  to  counterpropagate  through  a  medium,  causing  two 
holes  to  be  burnt  into  the  absorption  profile  due  to  the  interaction  with  opposite 
axial  velocity  components.  At  the  center  frequency  for  the  profile,  both  beams 
interact  with  the  same  velocity  subset. 24 
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At  all  other  points  there  are  two  separate  interaction  groups;  in  the  center,  the 
intensities  add  linearly,  and  there  is  only  one  point  of  interaction.  Because  of  this  single 
point  of  interaction,  there  is  a  dip  in  the  intensity  profile  of  the  beams  after  passing 
through  the  medium.  The  depression  in  the  intensity  curve  is  the  Lamb  Dip,  illustrated  in 
Figure  14  b.22 


Figure  14b.  The  Lamb  dip  in  the  intensity  profile.  The  dip  in  intensity  at  the  center 
is  created  because  the  beams,  which  previously  had  interacted  with  two  separate 
velocity  groups  in  the  medium,  now  influence  the  same  velocity  group.  The 
absorption  profile  rises  when  the  beams  cross  at  the  center  frequency  due  to  the 
mutual  interaction,  so  transmitted  intensity  dips.24 


The  interaction  of  counterpropagating  laser  beams  with  the  common  Zero  axial 
velocity  subgroup  is  the  key  to  the  understanding  of  Doppler-Free  Saturation 
Spectroscopy.  The  fact  that  there  is  common  interaction  for  both  beams  with  the  same 
molecule  velocity  component  allows  for  elimination  of  the  Doppler  information 
concerning  a  transition.  Once  Doppler  broadening  has  been  eliminated,  the  hyperfine 
interactions  can  be  easily  resolved. 
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D.  Doppler-Free  Saturation  Spectroscopy 


As  in  the  previous  section,  counterpropagating  two  beams  from  a  tunable  laser  source 
and  scanning  in  frequency  through  the  range  of  a  broadened  line  ensures  that  at  some 
point  both  beams  will  interact  with  the  same  group  of  molecules  in  velocity  space.  The 
general  diagram  used  to  design  Doppler-free  saturation  experiments  is  shown  in  Figure 
15.  Note  that  there  is  a  finite  angle  between  the  beams. 


Figure  15.  Typical  zero-Doppler  hole  burning  spectroscopy  setup.  The  two  beams 
counterpropagate  within  the  gas  cell,  and  cross  at  a  small  angle,  to  ensure  mutual 
interaction  with  the  molecules  at  the  center  of  the  absorption  profile.  Though  both 
beams  traverse  the  medium,  only  the  unchopped  beam  is  monitored  for  intensity 
fluctuation  caused  by  the  modulation  of  the  absorption  curve  by  the  presence  of  the 
modulated  beam.23 

There  are  vector  components  of  the  beams  that  are  counterpropagating  in  the  overlap 
region,  ensuring  the  mutual  interaction  with  the  axially  stationary  molecules.  The  weaker 
beam,  the  probe  beam,  is  the  beam  whose  absorption  will  be  measured  by  the  photo¬ 
detector  after  crossing  through  the  interaction  volume  in  the  cell.  The  stronger  pump 
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beam  is  modulated  at  a  reference  chopping  frequency  coth,  while  the  probe  beam  is 
allowed  to  propagate  as  a  continuous  wave  (CW).23 

The  pump  beam  intensity  oscillates  between  total  illumination  and  no  illumination, 
causing  the  absorption  profile  at  the  beam's  frequency  location  to  be  modulated.  A  hole 
appears  and  disappears  on  one  side  of  the  profile,  off  the  center  point.  At  the  same  time 
the  probe  beam  bums  a  smaller  hole  on  the  other  end  of  the  profile,  which  does  not 
oscillate.  Figure  16  shows  an  example  of  what  the  profile  may  look  like  for  a  given  off- 
center  frequency  component  in  the  absorption  profile  for  one  broadened  transition.24 


Figure  16.  Absorption  holes  for  0  Doppler  saturation  spectroscopy.  The  pump  hole 
oscillates  at  the  chopping  frequency,  while  the  probe  beam  burns  a  continuous  hole 
at  the  opposite  point  in  the  profile.  As  the  beams  scan  through  the  center  of  the 
profile,  the  CW  probe  beam  will  "feel"  the  suppression  of  the  absorption  curve 
caused  by  the  presence  of  the  pump  beam  and  will  have  a  transmitted  intensity  that 
is  modulated  at  the  chopping  frequency. 

As  the  laser  scans  in  frequency  space,  the  beams  eventually  will  interact  with  the 
same  velocity  subgroup  at  vz  =  0.  At  that  point  in  frequency  the  modulated  pump  hole 
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and  the  CW  probe  hole  overlap.  The  power  transmitted  by  the  probe  beam  will  be 
modulated  by  the  presence  of  the  chopped  pump  beam  in  the  same  spatial  location, 
interacting  with  the  same  velocity  group.  As  the  absorption  curve  is  suppressed  by  the 
pump  beam,  the  probe  beam  will  have  a  higher  transmission  through  the  medium.  On  the 
other  hand,  when  the  pump  beam  is  blocked  by  the  chopper  wheel,  the  probe  beam 
absorption  will  be  greater.  Figure  17  shows  the  modulation  of  the  pump  beam  by  the 
chopper  wheel,  with  the  corresponding  effective  modulation  on  the  probe  beam 
absorption  and  intensity. 
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Figure  17.  Probe  and  pump  beam  modulations 

A  lock-in  amplifier  can  be  used  to  detect  the  modulated  probe  beam  signal  created  by 
both  beams  being  present  in  the  same  spatial  region  at  the  same  time.  The  result  of  the 
chopping  of  the  pump  beam  and  the  subsequent  modulation  of  the  probe  beam  ensure 
spatial  resolution  for  the  Doppler-free  measurements.  The  lock-in  amplifier  turns  the 
oscillation  into  a  DC  signal  that  reflects  the  absorption  level  of  the  counterpropagating 
beams.  The  profile,  as  the  laser  frequency  is  scanned,  will  cross  through  the  zero 
Doppler  points  for  each  of  the  hyperfine  transitions,  creating  profiles  which  do  not 
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contain  any  of  the  Doppler-broadening  information.  The  lock-in  signal  maps  out  the 
naturally  broadened  hyperfine  structure  as  the  laser  scans  through  the  frequency  range  of 
the  broadened  line. 

E.  Doppler-Limited  Saturation  Spectroscopy 

There  is  another  way  to  ensure  that  the  two  beams  interact  with  the  same  velocity 
subgroup  within  the  velocity  distribution.  Instead  of  counterpropagating  to  cause  mutual 
interaction  with  the  vz  =  0  molecules  in  the  overlapping  spatial  region,  the  beams  can 
copropagate  to  overlap  not  at  the  Zero  Doppler  point  in  the  absorption  curve,  but  at  some 
other  velocity  component.  Since  the  beams  copropagate,  the  holes  are  burned  into  the 
profile  at  the  same  frequency  location,  but  away  from  the  Zero  Doppler  region  of  the 
profile.  The  interaction  between  the  pump  and  probe  beams  is  the  same;  the  difference  is 
that  the  lock-in  signal  now  maps  a  broadened  profile  as  the  laser  scans  in  frequency,  with 
all  Doppler  information  intact. 

The  spatial  overlap  of  the  beams  also  does  not  change  for  Doppler-limited  saturation 
spectroscopy.  The  probe  beam  will  still  only  show  a  measurable  signal  when  the  two 
beams  are  present  in  the  same  spatial  volume  at  the  same  time.  The  probe  beam  will  be 
modulated  at  the  pump  modulation  frequency,  as  before,  but  will  retain  the  complete 
Doppler  information  for  the  point  at  which  the  beams  overlap.  Scanning  in  frequency 
results  in  a  Doppler  profile  limited  to  the  volume  of  mutual  interaction  between  the  two 
copropagating  beams.  Figure  18  shows  the  conceptual  setup  image. 
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Figure  18.  Setup  diagram  for  Doppler-limited  hole  burning  spectroscopy.  This 
variant  of  the  standard  technique  allows  the  beams  to  copropagate  rather  than 
counterpropagate,  still  allowing  the  pump  beam  to  be  modulated  by  a  chopper. 


This  simple  variant  of  the  standard  technique  results  in  a  broadened  profile,  caused  by 
the  velocity  distribution  of  the  sample  gas  at  the  precise  spatial  location  where  the  beams 
overlap.  Therefore,  the  spatial  resolution  can  be  greatly  improved  over  typical  cell-length 
line  integrated  profiles,  where  the  beam  interacts  with  a  line  of  gas  throughout  the  whole 
length  of  the  cell.  Instead  of  losing  resolution  capability  because  of  the  Doppler  effect, 
the  hole  burning  technique  exploits  Doppler  information  to  allow  the  measurement  of 
either  the  naturally  broadened  hyperfine  component  lines,  or  the  spatially-resolved 
Doppler-broadened  profiles. 


HI.  Experimental 


A.  General  Setup 

The  primary  laser  source  for  the  experiment  was  the  Coherent  899  Ring  Dye  Laser, 
pumped  by  a  6  Watt  Argon  Ion  laser  at  the  5 14-nm  line.  The  ring  laser  was  setup  to 
operate  using  Rhodamine  590  dye,  with  an  output  beam  that  radiates  in  the  570-610  nm 
range.  This  range  determined  the  spectral  range  available  for  the  observation  of 
absorption  lines  in  molecular  iodine.  The  specific  transitions  chosen  were  well  in  the 
center  of  this  range,  to  allow  for  output  power  of  at  least  300  mW  from  the  dye  laser  in 
single  frequency. 

The  calibration  of  the  wavemeter  on  the  ring  dye  laser  assembly  was  very  important, 
as  at  times  systematic  alignment  problems  caused  offsets  of  almost  0.2  cm'1,  a  very 
significant  frequency  correction.  A  complete  calibration  procedure  is  offered  in 
Appendix  I,  by  using  Doppler- free  saturation  measurements  of  hyperfine  structure,  and 
comparing  measured  center  frequencies  for  the  naturally-broadened  lines  to  published 
values. 

The  ring  dye  laser  was  controlled  by  the  manufacturer's  autoscan  program,  which  was 
used  to  stack  scans  and  collect  data.  The  program  had  the  capability  to  run  scans  with 
various  wavelength  ranges,  resolution,  and  data  averaging  algorithms.  The  lock-in 
amplifiers  used  were  Stanford  Research  Systems  SR850s.  Chopping  frequencies  were  set 
to  3800  rotations  per  minute  (rpm's)  to  ensure  that  the  chopping  period  would  be  smaller 
than  the  set  lock-in  amplifier  time  constants. 

Two  separate  cell  designs  were  used  to  contain  the  iodine  gas:  the  first  was  an 
enclosed  glass  cell  with  one  inch  optical  windows  on  each  end,  in  which  the  true  pressure 
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was  unknown.  The  second  was  a  stainless  steel  housing  with  a  one-inch  optical  window 
fixed  to  one  end  and  three  ports  for  pressure  control/measurement,  and  gas  temperature 
determination  at  the  other  end.  The  experiment  was  actually  designed  in  two  phases, 
related  to  the  cell  design.  Phase  I  experiments  were  conducted  using  the  closed  glass 
iodine  cell,  which  was  primarily  utilized  to  obtain  a  rough  understanding  of  the 
measurement  technique,  and  to  prove  that  a  high  enough  signal  to  noise  ratio  (SNR) 
could  be  obtained  for  both  Doppler-free  and  -limited  saturation  measurements.  Phase  II 
experiments  involved  the  use  of  the  stainless  steel  cell,  in  order  to  better  control  the 
pressure  inside  the  cell  and  to  allow  for  comparison  to  reference  temperature 
measurements  made  by  a  thermocouple  within  the  same  cell  under  identical  conditions 

B.  Phase  I 

The  closed  I2  cell  was  first  used  to  demonstrate  that  the  technique  would  show 
enough  merit  to  justify  continued  investigation.  The  optical  table  was  first  setup  for 
Doppler-free  saturation  spectroscopy  measurements,  in  order  to  find  the  hyperfine 
splittings  within  given  transitions.  The  setup  is  shown  in  Figure  19.  For  this 
configuration,  a  pin  diode  (PD1)  measured  the  laser  power  to  ensure  that  fluctuations  in 
power  did  not  affect  the  absorption  profile.  The  photomultiplier  tube  (PMT)  measured 
the  side  fluorescence  from  the  iodine  cell,  so  that  the  broadened  line  could  be 
superimposed  with  the  hyperfine  structure.  PD2  was  used  to  measure  the  intensity  of  the 
probe  beam  after  it  passed  through  the  sample,  essentially  mapping  out  the  hyperfine 
spectrum  as  the  scan  progressed. 
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Figure  19.  Experimental  setup  for  Doppler-free  measurements  using  the  enclosed  cell 


One  important  thing  to  note  about  this  setup  is  the  fact  that  the  beams  were  aligned  to 
overlap  completely  throughout  the  iodine  cell.  The  two  adjustable  irises  ensured  near¬ 
perfect  overlap,  since  both  beams  were  adjusted  to  go  through  both  apertures,  set  at  the 
same  height  and  diameter.  This  allowed  for  very  strong  interaction  with  the  0  velocity 
gas  molecules  because  the  interaction  volume  was  large. 

The  PMT  measurements  allowed  for  large  scans  of  Doppler  broadened  lines  at  low 
resolution  to  be  made,  in  order  to  find  transitions  that  would  be  adequate  to 
representatively  demonstrate  the  technique.  The  first  two  lines  were  chosen  based  on 
their  isolation  in  frequency  space,  while  the  other  two  were  chosen  because  they 
consisted  of  overlapping  lines.  There  was  some  thought  at  the  inception  of  the 
experiment  that  overlapping  lines  might  provide  greater  temperature  sensitivity,  which 
justifies  the  exploration  of  transitions  that  overlap.  The  line  sets  chosen  were  those 
centered  at: 

(a)  17249.9623  cm1  (b)  17339.8187  cm1 

(c)  17340.2214  cm1  (d)  17340.7535  cm 1 

These  lines  had  lineshapes  that  are  shown  in  Figure  20.  They  are  presented  here  to 
illustrate  why  these  particular  transitions  were  chosen.  The  x-axis  on  these  diagrams  is 
the  frequency,  in  cm ‘and  the  y-axis  represents  relative  intensity. 
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Figure  20.  Broadened  Iineshapes  for  chosen  line  sets:  a)  17249.9623  cm'1, 
b)  17339.8187  cm1,  c)  17340.2214  cm1,  and  d)  17340.7535  cm1. 


The  identification  of  these  lines  was  determined  by  creating  energy  level  difference 
tables,  as  described  in  Chapter  II.  The  lines  were  attributed  to  the  following  transitions: 


Table  2.  Experimental  line  assignments 


v0  (cm1) 

Transition  (s) 

17249.9632 

P(66)  16-1,  P(51)  18-2 

17339.8187 

P(70)  17-1 

17340.2214 

P(53)  19-2 

17340.7535 

P(122)  16-0,  R(74)  17-1,  R(57)  19-2 

The  letter  designation  P  or  R  refers  to  the  rotational  branch  for  the  transition  in  Table 
2.  The  number  in  parenthesis  is  the  J"  number,  and  the  hyphenated  number  denotes 
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which  vibrational  level  change  occurred  in  conjunction  with  the  transition  between 
electronic  levels  (v'-v").  From  this  point  forward  in  this  document,  these  designations 
will  be  used  when  referring  to  specific  line  sets. 

Once  the  lines  were  chosen,  hyperfine  measurements  were  made  to  uncover  the 
Doppler-free  splittings  that  composed  the  broadened  lines.  A  scan  was  made  across  each 
of  the  transitions  with  the  following  scan  parameters: 

Scan  width  =  3  (GHz) 

Scan  interval  =  1  (MHz) 

Scan  rate=600  (seconds/10  GHz) 

The  hyperfine  spectra  were  recorded  for  the  four  transitions  and  the  setup  was 
changed  to  make  Doppler-limited  measurements  of  the  same  transitions.  The  two  beam 
components  were  redirected  to  copropagate  through  the  medium,  crossing  at  a  finite 
angle  0.  Figure  21  shows  the  Doppler-limited  setup.  The  mirror  (M2)  which  directed 
the  pump  beam  through  the  medium  was  mounted  on  a  translational  stage,  such  that  the 
lateral  position  at  which  the  beams  crossed  in  the  cell  could  be  changed,  while  leaving  the 
angle  the  same.  There  was  enough  freedom  to  move  the  beam  cross  point  from  one  end 
of  the  cell  to  the  other. 

To  create  a  heat  gradient  in  the  cell,  one  end  of  the  cell  was  wrapped  with  heater  tape. 
The  heater  tape  temperature  was  measured  with  a  thermocouple,  and  was  kept  constant 
after  equilibrium  was  reached  at  250  degrees  Celsius  (523  K).  With  this  cell,  it  was  not 
possible  to  have  thermocouple  temperature  measurements  of  the  gas  at  various  positions 
along  the  cell.  The  system  was  closed,  so  there  was  no  reference  temperature  distribution 
available  with  which  optical  results  could  have  been  compared. 
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Figure  21.  Experimental  setup  for  Doppler-limited  measurements  using  the  enclosed  cell 


Profiles  were  measured  at  the  1-inch  position,  at  the  3 -inch  position,  and  at  the  5 -inch 
position  in  the  6-inch  glass  cell.  Figure  22  illustrates  the  cell  length  and  the  positions 
where  the  beams  crossed  for  the  Doppler  limited  scans. 


Figure  22.  Cell  crossing  lengths  (shaded  region  is  heater  tape)  The  beams  as  shown 
are  copropagating  for  Doppler-limited  measurements. 

C.  Phase  II 

After  the  data  analysis  was  performed  on  the  data  from  Phase  I,  it  became  clear  that  a 
new  cell  design  was  needed.  It  was  imperative,  first  of  all,  to  be  able  to  control  and 
measure  the  pressure  in  the  cell,  to  avoid  the  possibility  of  pressure  broadening,  which 
could  have  added  to  the  linewidths,  creating  a  systematic  error  in  temperature 
measurements.  The  second  requirement  for  a  new  cell  was  to  have  the  ability  to  measure 
the  temperatures  at  positions  along  the  cell  in  order  to  have  a  comparison  standard  for  the 
Doppler  limited  optical  scans. 

For  convenience  the  cell  in  Figure  23  was  constructed  out  of  a  stainless  steel  housing, 
and  an  optical  window  was  affixed  with  epoxy  to  one  end  of  the  cell.  Three  ports  were 


available  at  the  other  end  of  the  cell.  The  two  side  ports  were  used  for  vacuum  system 
hookup,  and  for  the  placement  of  an  I2  gas  source  cell.  The  back  port  was  a  dual-purpose 
port  that  would  allow  interchange  of  either  a  second  window  or  a  thermocouple  fitting. 
The  ports  were  all  1/2"  in  outer  diameter,  so  that  they  could  be  fitted  with  vacuum-sealed 
Cajon  fittings.  The  thermocouple  fitting  for  the  end  port  allowed  a  probe  to  be  moved  to 
any  position  along  the  center  of  the  cell,  while  the  window  fitting  allowed  for  the  optical 
temperature  measurements. 


Window  Fitting 


Figure  23»  The  open  cell  design.  The  thermocouple  fitting  allows  a  thermocouple  to 
slide  through  the  center  of  the  cell,  to  measure  the  temperature  at  any  axial  point. 
The  window  fitting  replaces  the  thermocouple  fitting  at  the  end  of  the  cell,  to  allow 
tor  laser  propagation  and  optical  temperature  measurements. 
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Care  was  taken  in  this  phase  of  the  experiment  to  make  sure  that  the  hyperfine 
measurements,  the  Doppler-broadened  measurements,  and  the  reference  temperature 
measurements  were  conducted  under  as  close  to  the  same  conditions  as  was  possible. 

This  included  crossing  angles,  laser  intensity,  pressure,  and  relative  amplitude  ratios  of 
pump  intensity  to  probe  intensity. 

The  first  procedure  was  the  temperature  assessment  using  the  thermocouple. 

Pressure  was  kept  steady  at  0.2  Torr  +/-  0.002.  This  pressure  is  less  than  the  vapor 
pressure  of  I2  at  room  temperature  by  0.08  Torr,  so  it  is  lower  than  the  pressure  could 
have  been  in  the  closed  cell.  The  heater  tape  warmed  the  outside  of  the  cell  at  250 
degrees  while  the  cell  was  kept  at  vacuum.  The  gas  valve  was  opened,  left  at  pressure  for 
at  least  20  minutes  to  equilibrate,  and  measurements  were  taken  at  1,  2,  3,  and  4  inches 
along  the  cell  length.  The  schematic  for  the  thermocouple  measurements  is  in  Figure  24. 


Figure  24.  Reference  temperature  measurement  setup.  The  thermocouple  fitting 
was  used,  and  temperatures  were  measured  at  1, 2, 3,  and  4  inch  positions. 
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Instead  of  the  four  lines  investigated  in  Phase  I,  only  the  two  isolated  lines,  P(70)  17- 
1  and  P(53)  19-2  were  observed  in  Phase  II.  The  P(70)  17-1  transition  is  composed  of  15 
hyperfine  lines,  while  the  P(53)  19-2  transition  is  composed  of  21.  Sampling  these  two 
lines  accounted  for  all  the  possible  hyperfine  configurations  for  frequency-isolated  fines 
in  I2.  As  will  be  shown  in  the  next  chapter,  multiple  fine  sets  were  difficult  to  use  for 
precise  temperature  extraction,  although  there  is  still  evidence  that  suggests  overlapping 
lines  may  be  more  sensitive  to  small  temperature  variation. 

For  these  two  transitions,  hyperfine  measurements  were  made  in  the  same  manner  as 
in  Phase  I.  The  only  difference  in  the  two  setups,  as  is  evident  in  Figure  25,  is  that  with 
the  open  cell  the  beams  do  not  completely  overlap  at  a  zero  angle,  but  overlap  at  an  angle 
0  =  3.6  degrees.  This  angle  between  the  beams  was  necessary  to  keep  conditions  the 
same  for  Doppler  free  and  Doppler  limited  measurements;  if  the  beams  overlapped 
through  the  entire  cell  length,  then  the  Doppler  limited  spatially  resolved  measurements 
would  be  line-integrated  averages,  which  is  where  this  project  started. 

The  next  two  pages  show  the  setup  diagrams,  Figures  25  and  26,  for  both  the 
Doppler-free  and  Doppler-limited  experiments.  Doppler-limited  absorption  spectra  were 
obtained  by  crossing  the  beams  at  1,  2,  3,  and  4  inches,  in  order  to  compare  the  results 
with  the  thermocouple  measurements.  The  scan  parameters  and  cell  settings  were  the 
same  for  both  optical  procedures,  and  are  given  in  Table  3  below: 


Table  3  Scan  parameters  and  cell  settings  for  Phase  II 


Laser  Power  (prior  to  splitting) 

350  mW 

Scan  width 

5  GHz 

Heater  Tape  temperature 

250°  C 

Scan  interval 

1  MHz 

12  pressure 

0.2  Torr 

Scan  rate 

600  s/10  GHz 

Ratio  of  Probe/Pump  Intensity 

50% 
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Lens  ND:  Neutral  Density  Filter 

1:  Mirror  PD:  Pin  Diode 

:  Chopper  BS:  Beam  Splitter 
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Figure  25.  Doppler-free  configuration  for  open  cell 


Lens  ND:  Neural  Density  Filter 

1:  Mirror  PD:  Pin  Diode 

:  Chopper  BS:  Beam  Splitter 
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Figure  26.  Doppler-limited  configuration  for  open  cell 


D.  Data  Analysis  Procedure 


For  each  hyperfine  spectrum  measured,  the  data  was  transferred  into  ASCII  format 
and  imported  into  Peakfit  version  4.0.  Depending  on  the  J"  value  for  the  transition,  as 
explained  previously,  21  or  15  Lorentzian  lines  were  used  to  fit  the  data.  The  user- 
modified  Peakfit  functional  definition  and  settings  for  the  Lorentzian  hyperfine  lines 
were: 


ys=a3- 


aO 

It 

x-  (al  +-  a4) 

2 

a2 

where 

aO  =  line  amplitude  Shared 

al  =  line  center  frequency 

a2  =  natural  line  width  Shared 

a3  =  relative  line  amplitude  Locked 

a4  =  center  offset  frequency  Shared 


(19) 


Locked  values  indicate  that  the  user  will  define  values  for  each  peak;  these  values  do  not 
change  during  the  fitting  routine.  The  relative  amplitudes  were  locked  at  the  theoretical 
values,  proportional  to  2  F  +  1,  as  outlined  in  Chapter  II.  Shared  values  indicate  that  the 
fit  will  proceed,  ensuring  that  the  values  will  be  common  among  all  the  peaks.  The 
center  offset  shifts  all  the  lines  by  a  shared  amount  to  achieve  the  best  fit,  the  widths  are 
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necessarily  are  shared  value  because  the  lines  are  naturally  broadened,  and  the  overall 
amplitude  factor  simply  raises  or  lowers  all  the  peaks  to  match  the  relative  amplitude  of 
the  data  series. 

The  resulting  fit  then  gave  the  individual  locations  and  shared  FWHM  values  of  the 
hyperfine  lines.  The  Doppler  limited  data  for  the  same  transition  was  then  imported  into 
Peakfit,  and  the  information  from  the  hyperfine  scan  was  used  to  define  the  Doppler- 
broadened  fit  routine.  For  Phase  I  a  simple  Gaussian  profile  was  assumed.  The  user- 
modified  Gaussian  functional  and  parameter  settings  for  the  Gaussian  profile  on  Peakfit 
were: 


y=a3a0e 


1 

x—  (al-(-a4) 

2 

_  2 

a2 

where 

aO  =  line  amplitude 
al  =  line  center  frequency 
a2  =  Doppler  width 
a3  =  relative  amplitude 
a4  =  center  offset  frequency 


Shared 


Shared 


Shared 


Locked 


Locked 


(20) 
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For  Phase  II  a  Voigt  Profile  was  used,  so  that  the  Lorentzian  and  Gaussian  lineshapes 
could  be  individually  locked  or  shared,  as  necessary  for  the  fit.  The  modified  Peakfit 
tunctional  form  for  the  Voigt  Profile  and  its  parameter  setting  were: 


(21) 

where 


aO  =  line  amplitude 

Shared 

al  =  line  center  frequency 

Locked 

a2  =  Doppler  width 

Shared 

a3  =  natural  line  width 

Shared 

Locked 

a4  =  relative  amplitude 

Locked 

a5  =  center  offset  frequency 

Shared 

The  line  centers  and  the  Lorentzian  widths  were  locked  at  the  positions  resulting  from 
the  hyperfine  fit.  The  relative  amplitudes  remained  locked  at  2  F+l  proportions.  The 
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remaining  values  were  shared  among  the  peaks  for  the  fitting  routine.  The  amplitudes, 
again,  were  shared  to  account  for  changes  in  the  data  intensity  profile  and  the  center 
offset  compensated  for  uniform  shifts  in  the  frequency  axis  for  the  same  line.  The  shared 
Gaussian  width  and  its  fit  uncertainty  were  the  primary  factors  for  this  work,  since  the 
Doppler  effect  defines  the  width  of  the  Gaussian  component.  The  resulting  shared 
FWHM's  for  the  Gaussian  portions  of  the  lineshapes  were  used  to  extract  the  temperature 
information,  as  described  in  the  next  chapter.  For  all  cases,  the  fitting  routines 
considered  all  lines  simultaneously,  and  used  a  linear  baseline  subtraction. 

Signal  to  Noise  Ratios  (SNRs)  were  determined  for  both  Doppler-limited  and 
Doppler-free  data  by  taking  the  mean  peak  signal  level  divided  by  the  rms  deviation  of 
the  noise  around  the  mean  baseline.  For  this  reason,  when  SNRs  are  recorded  in  this 
document,  they  will  be  accompanied  by  the  mean  baseline  level,  as  well. 
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IV.  Results 


A.  Phase  I 

The  hyperfine  spectra  for  the  four  transition  sets-  P(70)  17-1,  P(53)  19-2,  the  triply 
overlapped  P(122)  16-0,  R(74)  17-1,  R(57)  19-2  states,  and  the  overlapped  P(66)  16-1, 
P(51)  18-2  states—  are  shown  in  Figures  27-30.  The  figures  show  the  component 
hyperfine  structures  behind  the  broadened  transitions  that  were  shown  in  Figure  20, 
measured  by  Doppler-free  saturation  spectroscopy  for  the  closed  glass  cell.  These 
spectra  are  representative  of  the  Phase  I  hyperfine  measurements,  which  were  repeated 
several  times  for  each  transition  under  varying  setup  conditions  to  try  to  obtain  a  good 
peak  SNR  for  each  transition.  With  the  almost  total  beam  overlap  in  this  configuration, 
the  interaction  volume  was  very  large  within  the  cell,  about  120  mm3  (120  x  1  x  1  mm) 
Therefore  decent  signaPnoise  levels  were  not  difficult  to  achieve,  and  are  recorded  in  the 
figure  captions. 

The  hyperfine  patterns  for  the  isolated  transitions,  P(70)  17-1  and  P(53)  19-2,  have 
patterns  of  15  and  21  hyperfine  lines,  respectively,  that  match  those  of  Figures  8  and  9. 
Since  these  are  the  only  two  possibilities  for  I2  transitions,  these  lines  are  representative 
transitions  for  the  purpose  of  this  experiment  and  will  provide  a  complete  data  set.  The 
transitions  that  overlap  with  each  other  may  in  the  end  provide  greater  sensitivity  to 
temperature  changes,  but  complications  to  be  discussed  later  arise  when  trying  to  fit  the 
relative  intensities  of  the  overlapped  lines. 
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Figure  27.  Hyperfine  spectrum  for  P(70)  17-1.  There  are  15  hyperfine  lines  that  comprise  this  transition.  Peak  SNR=18.68, 
mean  baseline  =  920.0  +/-  212.5. 


4000 
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mean  baseline  =  1986  +/- 


3000 


Figure  30.  Hyperfine  spectrum  for  P(66)  16-1,  P(51)  18-2.  There  are  2  overlapping  lines  that  comprise  this  set;  one  of  the 
lines  has  15  hyperfine  components,  and  one  has  21  hyperfine  components.  Peak  SNR=36.7,  mean  baseline  =  518  +/-  60.6. 


The  Doppler-limited  spectral  data  for  the  four  transitions,  measured  at  different  points 
along  the  non-uniformly  heated  cell  are  shown  on  the  following  pages.  For  these 
measuiements  the  Doppler-limited  setup  in  Figure  21  was  followed  using  the  closed  glass 
cell.  The  absorption  curves  for  the  three  different  temperature  points  are  superimposed 
on  the  same  plots  to  illustrate  the  differences  between  the  three  temperature  profiles.  It 
was  necessary  to  normalize  the  curves,  since  the  hot  end  of  the  cell  had  a  lower  number 
density  of  iodine  molecules  than  the  cold  end,  according  to  the  ideal  gas  law.  Pressure 
and  volume  in  the  closed  cell  were  fixed,  so  there  was  a  tradeoff  between  temperature 
and  number  density,  which  decreased  the  absorption  at  the  hot  end. 

Figures  31-34  are  presented  here  to  show  that  the  Doppler-limited  saturation 
spectroscopy  technique  does  satisfy  the  condition  for  spatial  resolution.  The  fact  that 
these  profiles  exist  demonstrates  that  the  absorption  profiles  are  spatially  resolved.  The 
CW  probe  beam  could  not  contribute  absorption  information  to  the  lock-in  signal  if  it 
were  not  interacting  with  the  same  spatial  volume  of  the  gas  as  the  modulated  pump 
beam. 

Line-integrated  measurements  within  the  same  cell  configuration  have  an  interaction 
volume  that  can  be  approximated  by  that  of  a  cylindrical  beam  with  a  1mm  radius 
interacting  through  the  entire  cell  length.  The  development  in  Appendix  B  provides  an 
approximation  for  the  interaction  volume  for  beams  crossing  at  an  angle  0=3.8  degrees. 
The  line  integrated  volume  was  approximately  Vu  =  120  mm3,  while  the  volume  for  the 
spatially-resolved  Doppler-limited  measurements  was  Vsr=12  mm3 . 
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Figure  31.  Absorption  profile  for  P(70)17-l  for  I2  at  3  Separate  Temperatures.  The  lightest  shade  indicates  the  lowest 
temperature;  increasing  darkness  means  increasing  temperature.  The  highest  and  widest  curve  is  the  hot  end 
measurement. 
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Figure  32.  Absorption  profile  for  P(53)19-2  for  I2  at  3  Separate  Temperatures.  The  lightest  shade  indicates  the  lowest 
temperature;  increasing  darkness  means  increasing  temperature.  The  highest  and  widest  curve  is  the  hot  end 
measurement. 


Figure  33.  Absorption  profile  for  P(122)  16-0,  R(74)  17-1,  R(57)  19-2  for  I2  at  3  Separate  Temperatures.  The  lightest  shade 
indicates  the  lowest  temperature;  increasing  darkness  means  increasing  temperature.  The  highest  and  widest  curve  is  the 
hot  end  measurement. 


Figure  34.  Absorption  profile  for  P(66)  16-1,  P(51)  18-2  for  I2  at  3  Separate  Temperatures.  The  lightest  shade  indicates 
the  lowest  temperature;  increasing  darkness  means  increasing  temperature.  The  highest  and  widest  curve  is  the  hot  end 
measurement. 


Though  the  main  objective  for  the  Phase  I  experiments  was  to  prove  the  spatial 
resolution  of  the  measured  absorption  transitions,  it  is  encouraging  to  note  the  dramatic 
differences  in  the  profiles  due  to  temperature  variation,  especially  for  the  overlapping 
lines.  The  isolated  lines  at  higher  temperatures,  as  expected  from  Figure  10,  broaden  into 
a  more  symmetrical  Gaussian-type  curve.  The  profiles  for  overlapping  transitions  tend  to 
flatten  out  for  the  hotter  temperatures,  giving  credence  to  the  notion  that  these  line  sets 
may  provide  better  temperature  resolution  than  the  isolates  fines. 

In  reference  to  the  Peakfit  functions  that  are  defined  in  Chapter  III,  Equations  19, 
20,  and  21,  the  relative  amplitudes  for  the  hyperfine  lines  are  set  and  locked  to  their 
theoretical  values  based  on  the  proportion  2F  +  1.  This  proportionality  factor,  however, 
only  relates  the  amplitudes  of  the  hyperfine  fines  within  the  same  transition.  For  the 
overlapping  transitions,  the  relative  heights  can  then  be  set  for  each  transition 
individually.  The  ratio  of  intensity  from  one  transition  to  another  is  a  completely 
different  story,  and  does  not  involve  a  simple  linear  relation. 

The  relative  amplitudes  for  two  different  transitions  are  temperature  dependent,  as 
hotter  temperatures  cause  the  increased  population  of  higher  ro-vibrational  states, 
changing  the  transition  rates23.  In  addition,  the  intensities  may  have  a  nonlinear  response 
to  incident  radiation  intensity,  such  that  the  amount  of  incident  light  energy  affects  each 
transition  differently.  This  would  mean  that  the  intensities  from  one  transition  to  the  next 
would  not  be  proportional  to  the  theoretical  intensities,  as  they  would  depend  upon  the 
laser  intensity.  Being  a  nonlinear  phenomenon,  the  prediction  of  overlapping  fine 
intensities  for  a  given  experimental  setup  would  be  a  nontrivial  task.  This  nonlinear 
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Relative  Intensity  (arb.  Units) 


response  could  possibly  be  temperature  dependent,  as  well,  further  complicating  the  issue 
of  accurate  temperature  extraction  using  hyperfine  structure  for  overlapping  lines. 

The  overlapping  transitions  therefore  comprise  a  very  complicated  system  for 
temperature  extraction  purposes.  Though  not  satisfying,  the  overlapping  transitions  are 
presented  in  the  previous  figures  merely  as  a  qualitative  illustration.  The  two  lines  which 
will  be  analyzed  for  the  quantitative  extraction  of  temperature  information  are  those 
which  are  isolated  in  frequency  space:  P(70)  17-1  and  P(53)  19-2  .  As  previously  stated, 
these  two  transitions  are  representative  of  all  possible  I2  transitions. 

The  hyperfine  fit  for  the  P(70)  17-1  transition  is  shown  in  Figure  35,  with  the 
corresponding  fit  parameters  given  in  Tables  4  and  5. 


Figure  35.  Hyperfine  fit  plot  for  P(70)  17-1.  This  is  the  same  spectrum  in  Figure 
27. 


65 


Table  4.  P(70)  17-1  hyperfine  fit  parameters 


Peak 

Center  (MHz) 

FWHM  (MHz) 

1 

-883.39 

9.94 

2 

-606.23 

3 

-591.76 

SNR 

4 

-587.78 

18.68 

5 

-573.63 

6 

-463.83 

r-squared 

7 

-451.97 

0.876 

8 

-440.24 

9 

-427.84 

10 

-299.75 

11 

-167.79 

12 

-163.80 

13 

-147.57 

14 

-145.28 

15 

-11.26 

The  resultant  value  for  the  hyperfine  linewidth,  9.94  MHz,  is  on  the  order  of  the 
expected  natural  linewidth.  A  signal  to  noise  ratio  of  18.64  is  adequate  to  differentiate 
the  hyperfine  lines  from  the  noise  floor.  The  r-squared  value  should  be  unity  for  an  ideal 
fit;  the  poor  value  comes  from  the  fact  that  the  addition  of  the  hyperfine  fines  do  not 
account  for  the  baseline.  In  the  figure,  it  is  also  evident  that  the  theoretical  intensities  do 
not  seem  to  account  for  the  two  hyperfine  lines  at  the  ends  of  the  spectrum;  this  will  be 
discussed  further  in  the  next  chapter. 

The  fine  positions  in  Table  4  were  used  to  set  the  centers  of  the  15  Gaussian 
components  to  fit  the  Doppler-limited  data.  Figures  36-38  show  the  Doppler-limited 
scans  and  their  accompanying  fit  curves.  The  hyperfine  Gaussian  curves  are  shown  at  the 
base  of  the  curves,  each  with  a  shared  Doppler  width  chosen  by  the  fit  routine.  The 
superposition  of  these  curves  accounts  for  the  solid  fine,  which  ideally  should  match  the 
trend  in  the  empirical  data  points. 
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Figure  36.  P(70)  17-1  Doppler-broadened  transition  at  the  cold  end  of  the  Phase  I 
cell,  measured  by  Doppler-limited  saturation  spectroscopy.  SNR=122.5,  mean 
baseline=18.5  +/- 13.6,  R2=0.993,  FWHM=581  MHz,  Extracted  Temperature=615  K 


Figure  37.  P(70)  17-1  Doppler-broadened  transition  at  the  warm  center  of  the  Phase 
I  cell,  measured  using  Doppler-limited  saturation  spectroscopy.  SNR=102.8,  mean 
baseline=24.0  +/- 13.8,  R2=0.986,  FWHM  =621  MHz,  Extracted  Temperature=707  K 
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Relative  Intensity  (arb.  Units) 


Figure  38.  P(70)  17-1  Doppler-broadened  transition  at  the  hot  end  of  the  Phase  I 
cell,  measured  using  Doppler-limited  saturation  spectroscopy.  SNR=37.5,  mean 
baseline=22.4  +/- 13.7,  R2=0.987,  FWHM=867  MHz,  Extracted  Temperature=1377 
K 


Table  5.  P(70)  17-1  data  summary  table 


W8BM 

m 

IB 

■iiiltill 

IMB 

Cold  End 

122.5 

0.993 

581 

615 

Warm  Center 

102.8 

0.916 

621 

707 

Hot  End 

37.5 

0.987 

867 

1377 

These  temperatures  are  much  higher  than  what  one  would  expect.  The  heater  tape 
itself  did  not  exceed  523  K,  as  measured  on  the  cell  surface,  so  the  lowest  optically- 
extracted  temperature  is  higher  than  the  heat  source.  Looking  at  Figures  36-38,  it  seems 
that  the  fit  curves  do  not  adequately  represent  the  data.  In  all  three  measurements,  the 
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parameterized  fit  tried  to  make  the  curve  wider  and  higher  than  the  data  suggests, 
resulting  in  higher  temperatures. 

Because  the  cause  of  the  obvious  error  in  the  optical  measurement  was  unknown,  the 
first  step  in  diagnosing  the  problem  was  to  make  sure  this  anomaly  was  not  transition 
dependent.  The  same  procedures  were  completed,  and  the  same  data  collected,  for  the 
P(53)  19-2  transition.  Figure  39  shows  the  Lorentzian  fit  to  the  hyperfine  spectrum  for 
this  transition,  and  Table  6  summarizes  the  main  fit  parameters. 


Figure  39.  Hyperfine  fit  plot  for  P(53)  19-2.  This  is  the  same  spectrum  shown  in 
Figure  28. 
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Table  6.  P(53)  19-2  Hyperfine  Fit  Parameters 


Peak 

Center  (MHz) 

FWHM  (MHz) 

1 

-604.25616 

12.28 

2 

-577.94301 

3 

-551.50746 

SNR 

4 

-327.7671 

69.6 

5 

-286.06186 

6 

-280.08122 

r-squared 

7 

-243.2564 

0.895 

8 

-161.52815 

9 

-143.76219 

10 

-131.69411 

11 

-119.66207 

12 

-23.136483 

13 

8.0834123 

14 

40.313806 

15 

118.02845 

16 

140.06365 

17 

160.64754 

18 

182.87031 

19 

281.94833 

20 

295.39243 

21 

310.31341 

The  FWHM,  though  slightly  higher  than  that  of  the  P(17)17-l  transition,  is  still  on  the 
same  order  of  magnitude  for  naturally  broadened  lines.  Again,  the  poor  r-aquared  value 
stems  from  the  fact  that  the  addition  of  the  hyperfine  lines  does  not  account  for  the 
baseline,  even  though  SNR  is  greater  than  in  the  last  case.  The  details  of  the  Doppler 
limited  analysis  for  P(53)  19-2  are  provided  in  Appendix  C.  Table  7  gives  the  summary 
of  the  data  for  this  transition. 


Table  7.  P(53)  19-2  data  summary  table 


ill! 

jBPj 

** 

HB81 

WKHKKk 

Cold  End 

53 

6.997 

422 

326 

Warm  Center 

51.2 

0.998 

441 

356 

Hot  End 

34.0 

0.989 

775 

1102 
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In  this  case,  the  warm  center  and  cold  end  measurements  are  more  reasonable,  but  the 
hot  end  temperature  is  still  a  nonphysical  result.  No  conclusions  can  be  drawn  at  this 
point  as  to  what  is  causing  the  discrepancy.  It  also  does  not  make  sense  at  this  point  to 
discuss  the  standard  deviation  calculated  from  the  fit.  The  Phase  II  measurements  were 
meant  to  eliminate  obvious  systematic  errors  in  order  to  narrow  the  possible  sources  for 
the  erroneous  temperature  measurement.  Phase  II  experiments  allowed: 

1)  Pressure  control  and  measurement  in  the  cell 

2)  Identical  conditions  for  the  Doppler-free  and  -limited  measurement  configurations 

3)  Knowledge  of  the  gas  temperature  at  the  measurement  points  along  the  cell 

4)  Use  of  a  Voigt  Profile  to  fit  the  broadened  lines 

In  essence,  the  systematic  errors  had  to  be  eliminated  to  ensure  that  the  discrepancy  in 
the  temperature  was  not  due  to  an  obvious  experimental  flaw. 
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B.  Phase  II 


The  first  data  collected  in  Phase  II  was  the  temperature  distribution  for  a  known 
pressure  of  gas  at  the  measurement  positions  in  the  new  cell.  The  thermocouple  was 
moved  to  each  new  position  and  was  allowed  to  equilibrate  at  pressure  for  twenty 
minutes  before  the  next  measurement  was  taken.  The  resulting  temperature  profile 
within  the  cell  is  shown  in  Figure  43,  and  is  presented  in  Table  8. 


Table  8.  Reference  temperature  values  at  each  measurement  point 


Temp(C) 

5 

33.3 

306.3 

4 

43.4 

316.4 

3 

57.8 

330.8 

2 

84.2 

357.2 

1 

130.3 

403.3 
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The  figures  and  tables  that  follow  summarize  the  data  analysis  for  Phase  II.  The 
analysis  was  the  same  as  for  Phase  I  measurements,  except  that  measurements  were  taken 
at  four  temperature  points.  The  hyperfine  spectra  were  measured  using  the  setup  diagram 
from  Figure  25  and  are  presented  below  in  Figures  41  and  42. 
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Figure  41.  P(70)  17-1  hyperfine  structure.  Doppler-firee  measurement  taken  with 
the  following  parameters:  ND=0.3  in  the  probe  beam;  Pressure=0.2  mm  Hg;  Lock- 
in  sensitivity=500  pV;  Lock-in  time  constant=30  ms;  SNR=55.24;  mean 
baseline=169.7  +/-  50.7 

This  is  a  larger  scan  in  frequency  than  was  completed  in  Phase  I.  The  longer  scan 
shows  that  the  hyperfine  lines  have  widths  that  add  to  produce  a  residual  broadened 
lineshape  at  the  base  of  the  hyperfine  scan.  This  trend  is  noticed  in  other  published 
documents  concerning  the  employment  of  Doppler-free  techniques  to  measure  hyperfine 
interactions.10,25  The  partitioning  of  data  for  Phase  I  partially  masked  this  trend,  which 
may  account  for  the  fact  that  the  natural  linewidths  that  were  fit  in  Phase  I  were  smaller 
than  those  in  Phase  II. 
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Figure  42.  P(53)  19-2  hyperfine  structure.  Doppler-free  measurement  taken  with 
the  following  parameters:  ND=0.3  in  the  probe  beam;  Pressure=0.2  mm  Hg;  Lock- 
in  sensitivity=100  pV;  Lock-in  time  constant=100  ms;  SNR=13.9;  mean 
baseline=471  +/-  179.9 


The  same  trend  is  noted  in  this  spectrum,  where  the  higher  energy  lines  all  tend  to 
have  base  widths  that  add  to  form  a  residual  broadened  lineshape.  It  is  not  clear  what 
could  be  causing  this  trend.  It  could  be  simply  that  the  natural  linewidths  are  broad 
enough  to  add  together  at  the  base,  or  that  some  other  broadening  mechanism  is  acting  on 
the  transition;  some  of  these  issues  will  be  addressed  in  the  next  chapter.  Figure  43 
shows  the  hyperfine  fit  for  P(70)  17-1,  with  relevant  parameters  listed  in  Table  9. 
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Figure  43.  Phase  II  hyperfine  fit  plot  for  P(70)  17-1.  This  is  the  same  spectrum 
shown  in  Figure  41. 


Table  9.  Phase  II  hyperfine  fit  parameters  for  P(70)17-l 


Peak 

Center  (MHz) 

FWHM  (MHz) 

1 

-251.2 

29.3 

2 

38.5 

3 

52.9 

SNR 

4 

56.9 

55.24 

5 

71.1 

6 

189.3 

r-squared 

7 

201.1 

0.937 

8 

212.9 

9 

225.3 

10 

353.3 

12 

494.5 

13 

498.4 

14 

514.7 

15 

517.0 
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The  resultant  hyperfine  linewidths,  29.3  MHz,  are  greater  than  those  found  in  Phase 
I,  as  stated  previously.  As  before,  the  two  end  lines  are  not  well-represented  by  the 
theoretical  intensities.  Appendix  C  contains  the  details  for  the  temperature  extraction,  the 
summary  of  which  is  presented  in  Table  10.  Figure  51  shows  the  comparison  between 
the  thermocouple  reference  measurements  with  the  optical  measurements  for  this 
transition. 


Table  10.  Data  summary  for  P(70)17-l  Phase  II  temperature  measurements. 


H 

: 

^■11 

m 

111— 

1 

IflMlII 

II 

4" 

143.8 

0.978 

774.3 

4.7 

1099 

316.4 

3  " 

139.3 

0.981 

792.6 

4.2 

1151 

330.8 

2" 

100.3 

0.981 

817.1 

4.2 

1224 

357.2 

1  " 

64.4 

0.98 

827.2 

3.4 
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Figure  44.  Temperature  comparison  for  P(70)17-l  transition.  Dotted  Line=optical 
measurements;  Solid  Line=thermocouple  measurements. 

Table  10  also  incorporates  the  statistical  standard  error  in  FWHM  from  the  fit,  simply 

to  indicate  that  it  is  small  when  compared  to  the  FWHM  value.  The  temperatures  for  this 
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transition  are  again  elevated,  when  measured  using  Doppler-limited  saturation 
spectroscopy.  Though  clearly  incorrect,  the  declining  trend  in  temperature  as  the 
measurement  position  changes  is  gradual  and  seemingly  continuous,  whereas  in  Phase  I 
there  were  disproportionate  jumps. 

Appendix  C  also  includes  a  study  of  the  fit  residuals,  which  are  remarkably  similar 
for  all  the  positions.  There  are  three  negative  peaks  and  two  positive  peaks,  at  about  the 
same  relative  position  in  frequency  space.  This  most  likely  points  to  incorrect  relative 
line  intensities,  since  there  is  no  way  to  attribute  a  sinusoidal  residual  variation  with  a 
baseline  subtraction  problem 

Figure  45  shows  the  hyperfine  fit  for  P(70)  17-1,  with  relevant  parameters  fisted  in 
Table  1 1 .  This  hyperfine  scan  is  a  closer  view  of  the  profile  in  Figure  42,  with  21 
Lorentzian  peaks. 
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Figure  45.  Phase  H  hyperfine  fit  plot  for  P(53)  19-2.  This  is  the  same  spectrum 
shown  in  Figure  42. 


Table  11.  Phase  II  hyperfine  fit  Parameters  for  P(53)19-2 


Peak 

Center  (MHz) 

FWHM  (MHz) 

1 

-610.8 

26.2 

2 

-584.5 

3 

-558.1 

SNR 

4 

-334.4 

13.9 

5 

-292.7 

6 

-286.7 

r-squared 

7 

-249.8 

The  26.2  MHz  natural  linewidths  are  again  greater  than  those  found  in  Phase  I,  but 
are  very  similar  for  both  transitions  in  Phase  II.  It  is  likely  that  the  true  natural  linewidths 
are  closer  to  30  MHz  than  it  is  to  10  MHz,  since  the  sectioning  of  data  for  Phase  I 
precluded  the  residual  broadened  lineshape  at  the  base  of  the  hyperfine  profiles.  Pressure 
broadening  was  not  influencing  hyperfine  lineshapes,  since  the  pressure  in  the  cell  was 
only  0.2  mm  Hg,  which  is  less  than  the  room  temperature  vapor  pressure  of  h,  which  is 
0.3 1  mm  Hg.  Generally,  pressure  broadening  does  not  show  significant  impact  until  a 
diatomic  gas  reaches  pressures  of  about  1  mm  Hg.23 

Appendix  C  contains  the  details  for  the  temperature  extraction,  the  summary  of  which 
is  presented  in  Table  12.  Figure  46  shows  the  comparison  between  the  thermocouple 
reference  measurements  with  the  optical  measurements  for  this  transition. 


Table  12  also  gives  the  statistical  standard  FWHM  error  from  the  fit;  again,  this  error 
is  far  too  minor  to  account  for  the  poor  temperature  determination.  Appendix  C  has  the 
accompanying  residual  study  for  all  four  measurements,  to  show  that  there  is  no 
overwhelmingly  obvious  trend  which  can  account  for  the  error.  There  are  many 
interesting  factors  about  these  results.  First,  the  more  careful  experimental  procedures  in 
Phase  II  did  not  change  the  erroneous  fit  values.  Though  a  systematic  error  is  still 
possible,  the  steps  in  Phase  II  should  have  eliminated  the  majority  of  the  unknowns  in 
this  area.  Second,  not  only  are  the  extracted  temperatures  wrong,  but  they  do  not  match 
one  another  between  the  two  transitions.  This  also  points  to  the  possibility  that  the  error 
may  not  be  systematic,  but  may  be  due  to  an  omission  in  the  theoretical  understanding  of 
this  technique. 

The  residual  trends  all  seem  to  take  on  the  same  character.  The  two-positive-peak, 
three-negative-peak  residual  after  the  fit  is  complete  is  consistent  for  both  transitions  for 
all  temperature  measurements.  This  in  itself  is  encouraging,  for  it  suggests  a  possible 
solution  to  the  temperature  error  problem;  however,  this  solution  is  not  clear  at  this  point. 
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V.  Discussion  and  Recommendations 


The  ability  of  Doppler- limited  saturation  spectroscopy  to  measure  absorption  spectra 
that  are  spatially  resolved  has  been  demonstrated.  The  nature  of  the  technique  ensures 
that  the  measured  interaction  occurs  only  at  the  spatial  location  where  both  beams  are 
present.  The  chopping  of  the  pump  beam,  and  its  presence  in  the  medium  modulate  the 
CW  probe  beam  to  produce  the  measurable  result.  At  the  points  where  the  beams  do  not 
overlap  the  probe  beam  is  not  modulated;  therefore,  the  resulting  lock-in  signal  is  zero. 
This  technique  offers  unique  capabilities,  as  it  can  be  used  in  any  Doppler-broadened 
medium.  The  spatial  resolution  of  such  measurements  is  a  marked  improvement  over 
line-integrated  optical  procedures,  allowing  for  measurements  in  a  1  x  1  x  12  mm  volume 
of  gas.  It  is  possible  to  improve  this  spatial  resolution  by  increasing  the  beam  crossing 
angle,  but  this  was  the  finest  resolution  that  could  be  achieved  within  the  physical 
limitations  of  the  sample  cell  that  was  used  in  the  present  experiment. 

The  overall  trends  and  lineshape  changes  are  encouraging,  for  they  indicate  that 
different  temperature  profiles  do  have  unique  characteristics.  Unfortunately,  the 
temperatures  extracted  from  these  spatially  resolved  profiles  did  not  return  realistic 
values.  There  appears  to  be  an  omission  of  a  detail,  probably  simple,  which  is  not 
allowing  the  procedure  as  designed  to  return  an  accurate  measurement  of  temperature. 
The  omission  can  be  either  in  the  experiment  or  in  the  theory. 

It  is  possible  that  the  discrepancy  is  still  due  to  some  systematic  error  with  the 
experiment.  Perhaps  the  number  of  parameters  allowed  to  float  in  the  Peakfit 
overextended  the  capabilities  of  the  package.  There  may  even  be  another  broadening 
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mechanism  which  was  unaccounted  for  by  the  experimenter  in  the  theoretical 
development.  Though  there  are  many  such  possibilities,  the  careful  design  of  the  Phase  II 
experiments  should  have  eliminated  the  majority  of  sources  that  could  have  caused  a 
400%  reading  error. 

Perhaps  there  is  an  omission  in  the  theory,  as  it  stands,  which  does  not  account  for  the 
temperature-dependent  broadening  of  the  fine  structure  lines  in  molecular  iodine.  The 
theoretical  intensities  that  are  related  to  the  factor  2F  +  1  did  not  really  match  the 
hyperfine  data  in  this  experiment,  especially  for  the  two  end  lines  in  P(53)  19-2.  The 
single  lines  at  either  end  of  the  spectrum  were  much  higher  than  some  of  the  doublet  or 
triplet  overlaps.  The  residual  broadened  lineshape,  too,  may  not  be  completely 
attributable  to  the  natural  linewidth.  This  residual  shape,  if  it  changes  with  temperature, 
could  drastically  affect  the  extracted  Gaussian  widths. 

The  foundationary  theoretical  development  for  hyperfine  interactions,  including  the 
intensity  proportionality,  was  formulated  for  microwave  frequencies.18  For  optical 
wavelengths,  it  is  possible  that  nonlinear,  frequency-dependent  effects  are  influencing 

line  intensities.  If  there  is  such  a  nonlinear  effect  occurring  with  relation  to  the  hyperfine 

/ 

intensities,  it  could  drastically  change  the  theoretical  understanding  that  is  the  foundation 
for  this  technique.  Unfortunately,  the  literature  that  exists  at  the  time  of  this  writing  only 
deals  with  the  positions  of  hyperfine  lines,  and  does  not  contain  empirical  data  which 
lead  to  an  understanding  of  the  intensities. 

In  addition,  if  saturation  effects,  which  were  not  considered,  are  occurring,  then  some 
line  intensities  could  be  suppressed  in  the  absorption  profile.  Such  effects  might  also 
explain  the  residual  lineshape,  due  to  the  high  absorption  within  the  medium,  and 


82 


possible  reemission  at  the  pump  beam's  chopped  frequency.  This  reemission  would  see  a 
medium  that  is  optically  transparent  because  it  is  saturated,  and  the  result  would  be  the 
Doppler-broadened  profile. 

Because  of  the  high  utility  of  this  spatially-resolved  temperature  measurement 
technique,  this  work  should  be  continued  in  order  to  narrow  down  the  field  of 
possibilities  for  producing  inaccurate  temperature  extraction.  A  prioritized  list  of 
possible  future  experiments  is  offered  in  Table  15. 


Table  13.  Future  Experimentation. 


— I  „  Ff  v  ^ 

Doppler-free  vs. 
Temperature 

If  hyperfine  widths  are  temperature  dependent,  then  this 
may  point  to  another  broadening  mechanism  affecting  the 
lines.  If  the  intensities  are  temperature  dependent,  then 
possible  theoretical  explanations  should  be  sought. 

Doppler-free  and 
Doppler-limited  vs. 
Crossing  Angle 

The  angle  at  which  the  beams  cross  affects  the  interaction 
volume,  and  so  may  point  to  a  saturation  phenomenon. 

Also,  off-angle  vector  elements  of  the  pump  beam's  wave 
propagation  vector  might  be  reaching  the  detector,  causing 
the  residual  lineshape. 

Doppler-free  and 
Doppler-limited  vs. 
Pressure 

Ensures  that  pressure  broadening  is  not  significant 

Doppler-free  and 
Doppler-limited  vs. 

Beam  Intensity  and 

Beam  Intensity  Ratio 

May  uncover  some  nonlinear  affects.  Also  may  show  that 
the  ratio  of  power  between  the  probe  and  pump  beams 
significantly  affects  the  resulting  spectra. 

Perform  Experiment  for 
other  Transitions 

May  uncover  some  nonlinear  affects.  Look  at  relative 
intensities  of  the  hyperfine  lines  and  compare  to  both 
theory  and  to  those  of  other  transitions. 

Perform  Symmetry 
Analysis  on  Doppler- 
Limited  Profiles 

May  increase  temperature  sensitivity  and  accuracy. 

Perform  Temperature 
Measurements  on 
Overlapping  Line  Sets 

Requires  reconciliation  of  relative  line  intensities,  but 
should  increase  temperature  sensitivity. 
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This  project  has  demonstrated  that  the  technique  shows  some  merit,  as  variations  in 
spatially-resolved  Doppler-profiles  have  been  detected.  Now  that  the  spatial  resolution 
and  the  temperature  differentiation  of  the  technique  have  been  demonstrated,  refinement 
of  procedures  and  theory  to  allow  an  accurate  extraction  of  temperature  will  complete  the 
picture.  There  is  not  as  of  yet  sufficient  data  to  offer  a  diagnostic  tool;  however,  this 
work  lays  the  foundation  for  future  endeavors. 
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APPENDIX  A 


Wavemeter  Calibration  Procedure 

This  procedure  is  meant  to  provide  guidance  for  making  very  fine  wavemeter 
calibrations  (+/-  a  few  MHz)  using  a  Ring  Dye  Laser  operating  in  the  570-610  nm  range. 
The  procedure  employs  Doppler-free  saturation  spectroscopy  to  measure  the  hyperfine 
spectra  for  several  transitions  in  molecular  iodine,  and  involves  comparison  of  specific 
lines  to  those  in  a  catalog  provided  by  Velchev  et.  al.25  Figure  A-l  shows  the  optical 
table  setup  diagram  that  accompanies  the  following  procedures: 


Figure  A-l.  Table  setup  for  ring  dye  laser  wavemeter  calibration 
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1.  Setup  the  optical  Table  as  shown  in  Figure  A-l,  using  a  glass  h  cell: 

a)  Ensure  the  beams  are  counterpropagating  and  completely  overlap  through  the  cell. 

b)  Ensure  LI  focuses  the  emitted  light  from  the  cell  to  the  PMT. 

c)  Set  the  chopper  frequency  very  high  (-3000  Hz). 

d)  Run  a  line  from  the  chopper  to  use  as  the  reference  for  two  lock-in  amplifiers. 

e)  Run  a  line  from  the  PD  to  use  as  the  signal  input  for  lock-in  #1 

f)  Run  a  line  from  the  PMT  to  use  as  the  signal  input  for  lock-in  #2 

g)  The  ring  dye  laser’s  output  should  be  at  least  350  mW  at  the  peak  of  the  dye  curve  in 
single  frequency  mode. 

Note:  lenses  may  or  may  not  be  necessary  to  focus  the  beams  through  the  cell. 
Overlapping  with  lenses  is  tricky  and  will  take  some  time.  Adjustable  apertures 
are  better  for  the  limiting  of  the  beam  diameter  going  through  the  cell. 

2.  Set  autoscan  for  a  100  GHz  Doppler-limited  emission  measurement.  The  results  of 
this  scan  will  be  compared  to  published  I2  reference  lines.  The  specific  center 
frequencies  are  given  in  the  atlas  by  Gerstenkom  and  Luc.5 

a)  Transfer  control  of  the  laser  to  the  autoscan  program. 

b)  Set  wavelength  =17334  wavenumbers. 

c)  Change  Data  Interval  to  50  MHz. 

d)  Change  Segment  Scan  Time  to  1 0s/ 1 0  GHz. 

e)  Change  Scan  Distance  to  100  GHz. 

f)  Edit  Sensor  List  to  only  read  from  lock-in  #2  (emission  measurements) 

g)  Execute  Scan 
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3.  Compare  resultant  spectrum  to  the  12  atlas.  The  lines  should  agree  with  some  degree 
of  reasonability  (at  least  0. 1  K.  If  not,  then  there  is  likely  a  misalignment  in  the 
wavemeter.  This  determines  a  rough  estimate  of  the  wavemeter  error. 

4.  Setup  autoscan  for  a  3.5  GHz  0-Doppler  hyperfine  absorption  measurement. 

a)  Set  wavelength  =one  of  the  start  scan  values  in  Table  A-l . 

b)  Change  Data  Interval  to  1  MHz. 

c)  Change  Segment  Scan  Time  to  500s/10  GHz. 

d)  Change  Scan  Distance  to  3.5GHz. 

e)  Edit  Sensor  List  to  only  read  from  lock-in  #1  and  2  (This  will  allow  you  to  display 
both  the  hyperfine  components  and  the  Doppler-broadened  line) 

f)  Execute  Scan 


Table  A-l.  Reference  positions  for  "t"  hyperfine  components25 

Transition  "t"  component  (Hz)  "t"  component  (K)  Start  Scan  (K)  Scan  Width  (GHz) 


P(65)  16-1 

5.171856956E+14 

1.725145786E+04 

17251.25 

3.5 

P(119)  18-1 

5.197223960E+14 

1.733607308E+04 

17335.85 

3.5 

P(1 16)  18-1 

5.199643760E+14 

1.734414466E+04 

17343.94 

3.5 

R(124)  18-1 

5.193047619E+14 

1.732214230E+04 

17321.93 

3.5 

5.  Compare  “t”  component  positions  to  obtain  a  refined  value  for  the  wavemeter  offset 
necessary  to  correct  the  error.  For  the  odd-J  lines  (P(65)  and  P(1 19)),  there  are  21 
hyperfine  lines,  while  for  the  even-J  lines  (P(l  16)  and  R(124)),  there  are  15  hyperfine 
lines.  For  even  J-lines,  the  "t"  component  is  defined  to  be  the  lowest  energy  line, 
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while  for  odd  J-lines,  it  is  the  second  lowest  energy  line.25  Figure  A-2  Clarifies  this 


definition. 


Figure  A-2.  Definition  of  the  "t"  component25 
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APPENDIX  B 


Estimating  the  Interaction  Volume  of  Two  Crossing  Cylindrical  Beams 

This  is  a  development  which  allows  for  the  approximation  of  the  volume  enclosed  by 
two  cylinders  intersecting  at  some  angle  0.  This  is  an  approximation  for  two  reasons:  1) 
TEM0,o  laser  beams  have  a  Gaussian  profile,  and  are  not  cylindrical;  2)  It  uses  simple 
geometrical  averages  to  provide  an  estimate  of  the  intersection.  As  a  first  approach,  we 
can  estimate  the  volume  inside  a  cylinder  by  averaging  the  volumes  of  two  rectangular 
beams.  One  of  the  rectangles  is  inscribed  within  the  cylinder,  while  the  other 
circumscribes  it,  as  in  Figure  B-l . 


Figure  B-L  Approximating  the  volume  of  a  cylindrical  beam  as  the  average  of  an 
inscribed  rectangular  beam  and  a  circumscribed  rectangular  beam.  The 
circumscribed  rectangle  has  a  cross  sectional  area  A2,  and  the  inscribed  rectangle 
has  a  cross  sectional  area  a2. 
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Knowing  the  radius  of  the  cylinder  (R),  we  know  the  lengths  of  the  sides  of  these 


rectangular  beams: 


Am  2-R 


(B.l) 


(B2) 


If  there  are  two  rectangles  crossing  at  an  angle,  the  area  of  intersection  is  a  polygon, 
which  is  shown  in  Figure  B-2,  with  sides  a,  p. 


Figure  B-2.  Polygon  created  by  the  intersection  of  two  rectangles.  A1  and  A2  are 
the  lengths  of  the  sides  of  the  square  cross  section  for  beams  1  and  2,  respectively. 


If  we  then  add  the  third  dimension,  the  volume  of  interaction  for  these  intersecting  beams 
is  then  the  area  of  this  polygon  multiplied  by  either  A1  or  A2,  whichever  is  smaller.  We 
now  have  a  way  to  describe  the  interaction  volume  of  two  independent  rectangular 
beams: 


V=  a-  P  •  sin(  0  )•  A^i 


(B3) 
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where 


a=A2/sin  (0)  P=  Al/sin  (0) 

To  approximate  the  volume  of  the  cylinder,  we  take  the  average  between  the  volume 
of  interaction  for  the  circumscribed  rectangles  and  the  inscribed  rectangles. 

Va=(V+v)/2  (B.4) 

where 

V,  v=the  interaction  volumes  for  the  circumscribed  rectangles  and  the  inscribed 
rectangles,  respectively 

This  approach  has  an  error  constant  associated  with  it,  which  is  simply  the  ratio  of  the 
approximate  volume  of  a  simple  cylinder  of  length  L  to  the  real  volume  of  the  same 
cylinder.  Any  volume  resulting  from  this  approximate  average  would  need  to  be 
multiplied  by  this  error  constant  in  order  to  provide  a  close  estimate. 

For  a  generic  cylinder,  radius  =  0.5  mm  and  length  =  10  mm,  the  approximate  volume 
is: 

A=1  mm  a=0.707  mm  V=10mm3  v=5  mm3 

Va=  15/2=7. 5  mm3 
And  the  actual  volume  is: 

Tir2L=7.854  mm3 

The  ratio  between  the  two  is  therefore  0.955.  This  ratio  is  true  regardless  of  which 
dimensions  one  is  to  pick.  So,  the  approximate  interaction  volume  must  be  multiplied  by 
this  factor  in  order  to  return  a  volume  close  to  that  of  two  intersecting  cylinders. 
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For  this  thesis,  the  two  beams  were  measured  at  equal  radii=0.5  mm.  The  beams 
crossed  at  an  angle=3.8  degrees.  From  these  parameters,  the  interaction  volume  is  found 
to  be  11.8  mm3. 
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APPENDIX  C 


Presentation  of  Additional  Data 

Phase  I.  P(53)  19-2  Temperature  Determination 

1.25  , - 


Av  (MHz) 

Figure  C-l.  P(53)  19-2  Doppler-broadened  transition  at  the  cold  end  of  the  Phase  I 
cell,  measured  using  Doppler-limited  saturation  spectroscopy.  SNR=53.0,  mean 
baseline=0.02  +/-  0.02,  R2=0.997,  FWHM  =422  MHz,  Extracted  Temperature=326  K 


Figure  C-2.  P(53)  19-2  Doppler-broadened  transition  at  the  warm  center  of  the 
Phase  I  cell,  measured  using  Doppler-limited  saturation  spectroscopy.  SNR=51.2, 
mean  baseline=13.7  +/- 12.7,  R2=0.998,  FWHM  =441  MHz,  Extracted 
Tempera  ture=356  K 
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Figure  C-3.  P(53)  19-2  Doppler-broadened  transition  at  the  cot  end  of  the  Phase  I 
cell,  measured  using  Doppler-limited  saturation  spectroscopy.  SNR=34.0,  mean 
baseline=14.4  +/-  9.7,  R -0.989,  FWHM  =775  MHz,  Extracted  Temperature=1102  K 

The  residuals  for  Phase  I  measurements  were  not  considered  in  order  to  focus  on  the 

elimination  of  systematic  errors  in  the  experimental  apparatus.  Phase  II  measurements 

will  include  a  fit  residual  analysis  to  determine  whether  or  not  there  might  be  a 

predictable  trend  in  the  data  error  between  the  measurements  and  the  fit  parameters. 
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Phase  II.  P(70)  17-1  Temperature  Determination  and  Residual  Study 


Figure  C-4.  P(70)  17-1  Doppler-broadened  transition  at  the  4  "  position  of  the 
Phase  II  cell,  measured  using  Doppler-limited  saturation  spectroscopy.  SNR=143.8, 
mean  baseline=0.0169  +/-  0.00696,  R2=0.978,  FWHM  =  774.3  MHz,  Extracted 
Temperature=1099  K 


Av  (MHz) 

Figure  C-5.  P(70)  17-1  residual  trend  for  the  4  "  Doppler-limited  fit. 


Figure  C-6.  P(70)  17-1  Doppler-broadened  transition  at  the  3  "  position  of  the 
Phase  II  cell,  measured  using  Doppler-limited  saturation  spectroscopy.  SNR=139.3, 
mean  baseline=0.208  +/-  0.00718,  R2=0.981,  FWHM  =  792.6  MHz,  Extracted 
Temperature=1151  K 


Figure  C-7.  P(70)  17-1  residual  trend  for  the  3  "  Doppler-limited  fit. 
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Figure  C-8.  P(70)  17-1  Doppler-broadened  transition  at  the  2  "  position  of  the 
Phase  n  cell,  measured  using  Doppler-limited  saturation  spectroscopy.  SNR=100.3, 
mean  baseline=+/-  0.00996,  R2=0.981,  FWHM  =  817.7  MHz,  Extracted 
Temperature=1224  K 


Figure  C-9.  P(70)  17-1  residual  trend  for  the  2  "  Doppler-limited  fit. 


Figure  C-10.  P(70)  17-1  Doppier-broadened  transition  at  the  1  "  position  of  the 
Phase  II  cell,  measured  using  Doppler-limited  saturation  spectroscopy.  SNR=64.4, 
mean  baseline=0.U4  +/-  0.0155,  R  =0.98,  FWHM  =  827.2  MHz,  Extracted 
Temperature=1254  K 


Av  (MHz) 


Figure  C-ll.  P(70)  17-1  residual  trend  for  the  1  "  Doppler-limited  fit. 
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Phase  II.  P(70)  17-1  Temperature  Determination  and  Residual  Study 


Figure  C-12.  P(53)  19-2  Doppler-broadened  transition  at  the  4  "  position  of  the 
Phase  II  cell,  measured  using  Doppler-limited  saturation  spectroscopy.  SNR=99.2, 
mean  baseline=0.567  +/-  0.101,  R  =0.99,  FWHM  =  555.8  MHz,  Extracted 
Temperature=566  K 
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Figure  C-13.  P(53)  19-2  residual  trend  for  the  4  "  Doppler-limited  fit. 
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Figure  C-14.  P(53)  19-2  Doppler-broadened  transition  at  the  3  "  position  of  the 
Phase  II  cell,  measured  using  Doppler-limited  saturation  spectroscopy.  SNR=68.0, 
mean  baseline=0.517  +/-  0.0147,  R2=0.99,  FWHM  =  643.4  MHz,  Extracted 
Temperature=759  K 


Av  (MHz) 


Figure  C-15.  P(53)  19-2  residual  trend  for  the  3  "  Doppler-limited  fit. 
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Figure  C-16.  P(53)  19-2  Doppler-broadened  transition  at  the  2  "  position  of  the 
Phase  II  cell,  measured  using  Doppler-limited  saturation  spectroscopy.  SNR=57.3, 
mean  baseline=0.205  +/-  0.0174,  R2=0.993,  FWHM  =  656.1  MHz,  Extracted 
Temperature=789  K 


Av  (MHz) 


Figure  C-17.  P(53)  19-2  residual  trend  for  the  2"  Doppler-limited  fit. 
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Figure  C-16.  P(53)  19-2  Doppler-broadened  transition  at  the  1  "  position  of  the 
Phase  II  cell,  measured  using  Doppler-limited  saturation  spectroscopy.  SNR=49.6, 
mean  baseline=0.0919  +/-  0.0202,  R2=0.995,  FWHM  =  672.2  MHz,  Extracted 
Temperature=828  K 


Av  (MHz) 

Figure  C-19.  P(53)  19-2  residual  trend  for  the  1"  Doppler-limited  fit. 
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